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Abstract. We present theorems which provide the existence of invariant 
whiskered tori in finite-dimensional exact symplectic maps and flows. The 
method is based on the study of a functional equation expressing that there is 
an invariant torus. 

We show that, given an approximate solution of the invariance equation 
which satisfies some non-degeneracy conditions, there is a true solution nearby. 
We call this an a posteriori approach. 

The proof of the main theorems is based on an iterative method to solve 
the functional equation. 

The theorems do not assume that the system is close to integrable nor 
that it is written in action-angle variables (hence we can deal in a unified way 
with primary and secondary tori). It also does not assume that the hyperbolic 
bundles are trivial and much less that the hyperbolic motion can be reduced 
to constant. 

The a posteriori formulation allows us to justify approximate solutions pro- 
duced by many non-rigorous methods (e.g. formal series expansions, numerical 
methods). The iterative method is not based on transformation theory, but 
rather on succesive corrections. This makes it possible to adapt the method 
almost verbatim to several infinite-dimensional situations, which we will dis- 
cuss in a forthcoming paper. We also note that the method leads to fast and 
efficient algorithms. We plan to develop these improvements in forthcoming 
papers. 

Keywords: whiskered tori, hamiltonian systems, small divisors, KAM theory 
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1. Introduction 

The goal of this paper is to prove some results on persistence of invariant tori 
for symplectic and exact symplectic maps and flows. 

We will assume that the motion on the torus is a Diophantine rotation and that 
the remaining directions are as hyperbolic as allowed by the symplectic structure 
(if the remaining directions are not void such tori are commonly called whiskered 
tori). 

More precisely, as it is well-known, the preservation of the symplectic structure, 
together with the fact that the motion on the torus is a rotation, implies that the 
symplectic conjugate direction to the tangent of the torus is not hyperbolic. We 
will assume that the remaining directions in the tangent bundle of the phase space 
at the torus are spanned by a basis of vectors which contract exponentially in the 
future or in the past. 

To make the previous statements more precise, we discuss first the case of maps. 
As we will show, results for flows can be readily deduced from the ones for maps. 
Given an exact symplectic map F from an exact symplectic manifold (M , O = da) 
into itself (for the purposes of this preliminary exposition, we will take A4 to be an 
Euclidean manifold, even if we will indicate how to eliminate this restriction later), 
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and a frequency vector G R , we seek an embedding K -.r ^ M. satisfying 

(1) {FoK){e)^K{9 + Lu), eeT^^R^/zK 

Equation ([T]) implies that the range of K is invariant under F. If ii' is an 
embedding, we obtain that K{T^) is a torus contained in A4, invariant by F and 
that the dynamics on it is, up to a change of coordinates, just a rotation of rotation 
vector Lu. 

The main resuh of this paper will show that if we can find a function K which sat- 
isfies some non-degeneracy assumptions and which satisfies ([T]) up to a sufficiently 
small error, then there is a true solution nearby. 

Differentiating the functional equation ([J) with respect to one gets 

DF{K{0))DK{e) = DK{e + uj). 

Geometrically, this shows that the tangent vector-field DK{9) is invariant and does 
not grow or contract under iteration of the action by the map. 

As we will see in more detail in Section 1423] if the map preserves the symplectic 
form Q and K is a solution of the invariance equation, there exists an analytic matrix 
valued function A{9), such that 

(2) DF{K{e))[J{K)-^DKN]{e) = DK{9 + lu)A{0) + [J{K)-^ DK N]{9 + u)), 
where J is the matrix representation of the symplectic form and 

N{e) = [DK{9)'^DK{d)]-\ 

As a consequence, [J{K)^^ DK N] [6) cannot grow more that polynomially. Hence 
we obtain that the center subspace of Tx(9)M is at least a 2?-dimensional space 
spanned by DK{e) and [J{K)-^DK N]{9) (we wifi show that range DK{9)n range 
[J{K)~^DK N]{9) = {0} because the image of the torus is a co-isotropic manifold.) 

For approximately invariant systems, the previous identities are just approximate 
and this implies that the center direction is at least 21. We will assume that indeed 
the dimension of the center subspace is exactly 21. That is, we will assume that the 
tori we consider are as hyperbolic as allowed by the fact that the motion on them 
are rotations and that the system preserves the symplectic structure. 

The main non-degeneracy assumptions on the approximate solution are a) that 
the other directions in Tx(e)M are hyperbolic. That is, they are spanned by vectors 
which contract exponentially fast in the future or in the past, b) that there is some 
twist condition, that is, that the matrix A in ([2]) is invertible. 

We will use a KAM iterative method to show that, if we are given a function K 
which solves ([T]) up to an error which is sufficiently small with respect to the prop- 
erties of the non-degeneracy conditions a), b) above, then there is a true solution 
close to this approximate solution. 

These results based on validating an approximate solution — which we call a 
posteriori — imply the usual persistence results (one can take as approximate solu- 
tion of the modified system the exact solution for the original one). Nevertheless, 
the a posteriori results can be used for other purposes. For example a posteriori 
results can be used to validate solutions obtained through any method such as nu- 
merical approximations or asymptotic methods. The validation of Lindstedt series 
leads to estimates on their domain of analyticity. The paper |Mas05| for instance 
considers Lindstedt series of whiskered tori. 
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A posteriori results also lead automatically to Lipschitz dependence on para- 
niaters and, with a bit more of work, to difFerentiable dependence on parameters. 
The a posteriori approach to KAM theorem was emphasized in JMos66b| IMos66a[ 
IZeh751 rZeh76] . There, it was pointed out that this a posteriori approach automati- 
cally allows to deduce results for finitely differentiable systems. We refer the reader 
to [dlLOlc] for a comparison of different KAM methods. 

In the present paper we deal with finite-dimensional maps and flows. In the 
forthcoming second part of it we consider coupled map lattices [FdlLSOS] . The case 
of Partial Differential Equations, which can be treated in a similar way but involves 
technical difficulties, is postponed to a forthcoming paper [dlLSOT] . 




Results on whiskered tori similar to the finite-dimensional ones of this paper have 
been considered several times in the literature. The first ones are |Gra74[ FZehTGj . 

The approach in [Zeh76j — which also takes the a posteriori format — is based on 
[Zeh75| which consists of finding a change of variables which reduces the system to a 
normal form which obviously possesses an invariant torus. This change of variables 
is accomplished by applying a sequence of canonical transformations. The method 
of proof introduced here is not based on successive transformations but rather 
on successive corrections introduced additively. This makes the estimates easier 
to establish and it leads to efficient numerical implementations. In order to be 
able to solve the equations, we take advantage of some cancellations due to the 
preservation of the symplectic structure that were also pointed out in |JdlLZ99l 
HLG JV05, IdlLOlcl 

The method of |Zeh76| proves the result for periodic Hamiltonian flows. The 
result for diffeomorphisms is proved in |Zeh76| by interpolating diffeomorphisms by 
periodic flows and then applying the results for periodic flows. The proof we present 
here proceeds along the opposite route. We prove first the result for diffeomorphisms 
and, then, deduce the result for flows taking time-one maps. Giving a direct proof of 
the result for persistence of whiskered tori for maps has been suggested as somewhat 
desirable in J. Moser's Mathematical review for [Zeh76j . We also provide such a 
direct proof. Of course, if one uses normal forms — as in [Zeh76| — it is natural 
to consider flows since the normal forms require only the study of the Hamiltonian 
function, which transforms very well. In the method presented here, the geometric 
cancellations are much more transparent in the case of diffeomorphisms. 

Among other results for finite-dimensional systems, we call attention to [ValOO| . 
which uses a method similar to that of |Arn63| . The paper [ValOO] has the advantage 
that it is a first order method (i.e. that each step of the Newton iteration requires 
to solve only one small divisors equation). As a consequence, the size of the gaps 
among tori in near integrable systems, the loss of regularity as a function of the 
Diophantine exponent and the required minimum regularity are smaller than these 
of the second order methods. A comparison between first and second order methods 
to prove KAM results can be found in [dlLOlc] . The paper |Sor02| (see also the 
sketch in |Moe96j ) uses a reduction to a normally hyperbolic manifold and then 
applies the standard KAM theorem for Lagrangian tori. Of course, since normally 
hyperbolic manifolds are in general only C' , the above method cannot produce C°° 
or analytic tori. On the other hand, we note that the method of |Sor02| IMoe96| 
leads to very good regularity conclusions for finite differentiable systems and also 
to good estimates on the measure occupied by the tori. We also call attention 



to [ZLLOSl IHLY061 ILYOSbi ILYOSal ISev06l ISevQQl lEIiOTl IeIi89] which consider also 
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tori with hyperbolic and eUiptic directions and relax the twist conditions and the 
differentiabihty requirements. The paper [GG02j considers analytic perturbations 
which depend only on the angles of reducible tori satisfying a twist condition and 
uses a direct resummation method. 

As compared with previous finite-dimensional results, the method presented here 
has the advantage that one does not need to assume that the hyperbolic bundles 
are trivial (and much less that the motion in the hyperbolic directions is reducible 
to a constant linear map). Tori with non-trivial invariant bundles appear naturally 
in one parameter families after crossing a resonance, see |HdlL07j . 

Also, we do not need to assume that the system is given in action-angle coordi- 
nates, something which is convenient if we are working in situations when the action- 
angle coordinates are singular. For instance, in the study of diffusion one is lead 
naturally to the study of whiskered tori near resonances (see [DdlLSOS, DdlLSOGj ). 
In this case, the action-angle variables are singular and avoiding its use leads to 
better estimates. 

For symplectic ODE's we will also prove a translated tori theorem. From this 
general version we will deduce the results for exact symplectic ODE's using a van- 
ishing lemma. We note that the approach of proving a translated torus theorem 
was introduced in [Rus76a in the one degree of freedom case. 

The method presented here lends itself to a very efficient numerical implemen- 
tation (see [HdlLSOS] ). The only functions to be considered are functions with a 
number of variables equal to the dimension of the torus itself (independently of 
the number of variables of the ambient space). Of course, when studying infinite- 
dimensional systems — PDE's or coupled map lattices or chains of oscillators — 
studying functions with the number of variables of the phase space is prohibitive. 
When implementing our method, if we discretize the tori by N Fourier coefhcients, 
the algorithm presented here only requires storage of order of N and a Newton step 
takes only order of N log{N) operations using the Fast Fourier Transform. This 
seems to be significantly faster than other algorithms. Actual implementations are 
now being pursued and will be the subject of a forthcoming paper (see [HdlLSOS] ). 
We refer the reader to [HdlL06b. .HdlL06ai, IHdlLOT] for analysis and implementation 
of related algorithms. 

2. Definitions and notations 

Before presenting the basic ideas and the results of our method, we introduce 
some notations and definitions which are useful for our purposes. All definitions 
are rather standard and we collect them here mainly to set the notation. 

2.1. Diophantine vectors. In the study of invariant tori one needs an arithmetic 
condition over the frequency vector. In the case of maps the notion of Diophantine 
vector is the following. 

Definition 2.1. Given k > and v >l, we define D{k,v) as the set of frequency 
vectors G M' satisfying the Diophantine condition: 

\uj ■ k - n\'^ < Klkl" , for alike I.^\{0} and n £ Zi, 

where ■ means scalar product, \k\ — \ki \ + ■ ■ ■ + and ki are the coordinates of k. 
We will say that to £ D{k, v) is Diophantine. 

For vector-fields the corresponding notion is the following. 
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Definition 2.2. Given k > and v > I — 1, we define Dh{n,i^) as the set of 
frequency vectors cj G M' satisfying the Diophantine condition: 

|w • fc|-i < Klfcl", for all k e \ {0} 

with the same notation as in Definition \2.1\ 

The two conditions are closely related since lo G Dh{K, v) with oji ^ if and only 
if {u>2/i^i, ■ ■ ■ uji/uJi) £ D{k\ u) for some k' . The geometric and measure properties 
of the sets of Diophantine vectors have been extensively studied. These results 
translate immediately into statements about the abundance of KAM tori. 

Definition 2.3. Let T' = R'/Z' and f G ^^(T'). We denote avg (/) its average 
on the l-dimensional torus, i.e. 

avg(/) = / f{e)de. 

Jv 

Definition 2.4. Given G we introduce the rotation over T' of rotation vector 
ui: 

TU9) = e + Lj. 

2.2. Functional spaces, functions and operators. We will denote Dp the com- 
plex extension of the torus of width p, i.e. 

(3) Dp = {ze C'/Z' I llmz^l <p, i = l,...,l}. 

We denote by | • | the supremum norm on or C^. The sup norm makes several 
estimates independent of the dimension of the manifold, which are useful when 
considering infinite-dimensional problems. However on Z' we will use the norm 
fc| — |fci |+...+|fc; |. Furthermore, for finite differentiability purposes, we consider the 
following norms: given g analytic, with bounded derivatives in a complex domain 
B, and m G N we introduce the following C™-norm for g 

Ifflc^CB) = sup sup|D''5(z)|. 

0<|fc|<m zeB 

Let Ap be the set of continuous functions on Dp, analytic in the interior of Dp 
with values on a manifold A4, which is assumed to be Euclidean. We endow the 
space Ap with the usual supremum norm 

= sup \u{z)\. 

z£Dp 

We have that {Ap, \\ ■ \\p) is a Banach space. In particular, ||u||o = 1|u1|2,oo(ti). 
We also recall the following convexity property (see |Rud74i Lemma 12.8]). 

Proposition 2.5. Let < pi < p2 and assume that f G Ap^. Then, for every 
G [0, 1] we have: 

(4) ll/ll.p.+(i-.)p. < ll/ll^jl/ll^;". 

In particular, taking pi = and 6 = 1/2, 

(5) ||./t./2<||/||^/i(TMll/llp^'- 
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We will also consider spaces of continuous functions on Dp analytic in its interior 
and taking values on finite-dimensional vector spaces, for instance in spaces of 
matrices. When endowed with the supremum norm, these function spaces are also 
Banach spaces. 

In particular, we will also need some norms of linear maps on the tangent space 
TK{e)M. with K{6) e AA, where K is an embedding. More concretely, let A{9) 
be a continuous linear operator from Tfc(g)M. into itself depending on the variable 
9 e Dp. Then we define ||A||p by 

||A||p = sup sup \\Ai9)v\\p. 

3. Setting of the problem and results 

3.1. Geometric setup. We will consider the Euclidean manifolds A4 = M.'^'^ and 
— R2d-2i X M' X T'. In the second case we can consider the universal covering 
R^'' of M and lift the maps defined on M to maps F, defined on M?"^, such that 
ttF = Ftt, where tt : M^"* — > A4 is the canonical projection. Even if we pass to 
the covering we will use the symbol A4 to refer to the manifold. These manifolds 
obviously admit complex extensions by considering M C C and T = R/Z C C/Z. As 
we will see, these different possibilities are convenient when we consider tori whose 
embeddings are topologically different. For example, tori which are contractible to 
tori with different dimensions. We will use the same symbol A4 for the complex 
extension of the manifold or its covering. 

For convenience of notation, we will endow these manifolds with the standard 
Riemannian metric, even if this may not be natural for the problem at hand. For 
us, the metric will only play the role to measure sizes and therefore any equivalent 
metric will give a similar result. The standard metric will have the advantage that 
it will allow us to use matrix notation for adjoints. In matrix notation, thinking 
of vectors as column vectors, we can write a^b = {a,b). On the other hand, we 
note that the length of vectors will always be the supremum norm and the norm of 
matrices will be the operator norm associated to the supremum norm on vectors. 
Of course, for finite-dimensional problems the supremum norm is equivalent to the 
Euclidean norm. 

We will assume that the Euclidean manifold A4 has an analytic exact symplectic 
form Q with primitive a, i.e. Q = da. For each z G M,\et J{z) : T^M — > T^M be 
the isomorphism such that 

mv) ^<^,j{z)ii>, 

where <, > is the Euclidean product on T^AA. 

We will not assume that J{z) has the standard form. We do not assume either 
that J induces an almost-complex structure on TAi. This generality is useful in 
some applications, (celestial mechanics, numerics, ...) when we use some system 
of coordinates — e.g. polar coordinates — which lead to non-standard symplectic 
matrices. 

Remark 3.1. As we will see in the proof, we are not using much the Euclidean 
structure of the manifolds. In Section FTBl we will present the modifications needed 
to work on other manifolds. 

More precisely, we will show that it is possible to work out the proof in a neigh- 
borhood U of the zero section of a bundle E" ® E" ® E"^. The bundle E" wiU be 
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shown to be trivial (as a consequence of the fact that the motion on the torus is a 
rotation, the preservation of the symplectic structure and the fact that the dimen- 
sion of the center space is 21, see Section [4.2.21) . but the others bundles — which 
correspond to the hyperbolic directions need not be trivial. 

We note that equation ([T]) is geometrically natural since it can be formulated in 
any manifold. 

In the following write up the Euclidean structure enters in two ways: one is 
a purely notational one and can be eliminated at the price of a typographical 
nightmare. When we only have approximate solutions, we will denote the error 
just as FoK{6) — K{9-\-lo) rather than eyi]y^^g^^^{F o K {9)) . We will also compare 
vectors in TpaK{9)M. with vectors in TV(g_|_^)A^. This can be done by introducing 
connectors as in |HPPS70j . so that what we denote DF o K{6 + uj)DF o K{9) is 
really S^^°^fll'^^^ DF o K{9 + uj)S^°f_^l\DF o K{9). See Definition EH and the 
discussion after (particularly equation (|17p ). 

A second, and more serious way in that the Euclidean space enters is that, 
to implement the iterative step in KAM theory, we will use Fourier series. This 
certainly requires that the functions take values in a vector space. Fortunately, 
this happens only in the center directions. In the hyperbolic directions there are 
geometrically natural ways to solve the iterative equation. This is why we are 
requiring that the center bundle is trivial, but we do not need the triviality of the 
hyperbolic bundles. 

Of course, the fact that we work in a set U as above is no loss of generality 
because, if there is a whiskered torus, by the tubular neighborhood theorem, we 
can identify a neighborhood of the torus with a neighborhood of the zero section 
of the normal bundle. 

3.2. Setting of the problem and results for maps. The main purpose of the 
theory we are going to develop is to construct invariant tori for exact symplectic 
maps. We recall the following 

Definition 3.2. Let {A4, il = da) be an exact symplectic manifold. A map F from 
M. into itself is exact symplectic if there exists a smooth function W on M such 
that 

F*a = a + dW. 

In particular, every exact symplectic map is symplectic, i.e. F*fl — Q. 

Heuristically, our problem is the following: let F be an exact symplectic map 
and uj G D{k,,u). We want to construct an invariant torus for F such that the 
dynamics of F on it is conjugated to T^. To this end, we search for an embedding 
K : Dp D T' ^ in Ap such that for all 9 G Dp, K satisfies the functional 
equation 

(6) F{K{9)) = K{T^{9)). 

Notice that if ^ is satisfied, the image under F of a point in the range of K will 
be also in the same range. Hence, since K is an embedding, the range of K will be 
an invariant torus. 

The assumptions of our results will be that we are given a mapping K that 
satisfies ([6]) up to a very small error and which satisfies some non-degeneracy and 
hyperbolicity assumptions. We will prove that then, there is a true solution of ([6]) 
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close to K. Wc will also prove that the solution of ([6|) is unique up to composition 
on the right with translations. 

The exactness of the map F is important for the existence of a solution to © . 
It is easy to construct examples of symplectic non-exact symplectic maps without 
invariant tori. For instance, consider = T x M with the standard symplectic 
structure. The translation in the M-direction is a symplectic non-exact symplectic 
map without any invariant torus. 

To construct the desired invariant torus, we consider a parameter A G M' and 
introduce a translation term in equation ^ depending on 9. 

We then consider the following functional equation, where G is a suitably chosen 
function of 9 taking values in 2d x I matrices and whose unknowns are both K and 
A 

(7) F{K{9)) + G{9)X = K{TU0))- 

The introduction of this parameter A will allow us to sidestep several technical 
complications and then we will show that, since F is exact symplectic, the geometry 
implies that A = 0. The fact that the dimension of the parameter A is / is important 
for our purpose. We also mention that it is possible to use the parameter A to 
weaken non-degeneracy conditions by taking A G R^' instead of R'. In such a case, 
G is a. 2d X 21 matrix. 

Remark 3.3. The introduction of the parameter A is also motivated by numerical 
calculations (see [HdlLSOS] ) . It leads to more stable computations. More impor- 
tantly, it is useful in the numerical computation of secondary tori (i.e. tori generated 
by resonances, which have some contractible directions). 

We go through a KAM technique to prove the existence of such a pair {X,K). 
To this end, we introduce the operator J-'^^ 

(8) J^^{X,K)^FoK + GX-KoT^, 
where 

(9) G^[J{Ko)-'DKo]oT^ 

is a function defined on T' and where Kq stands for an approximate whiskered torus. 
We will write G instead of its explicit form in many of the following results. As we 
will see later, the important property of G is that translations along the direction 
of G can change the cohomology of the pushforward in the center directions. 

The method is based on a careful study of the linearization (around a given pair 
(A, K)) of the operator T^j. We will show that this linear operator is approximately 
invertible in a suitable sense. 

For that, we have to introduce several non-degeneracy conditions. 

Definition 3.4. Given A G and an embedding K : Dp D T' ^ 7W we say that 
the pair {X,K) is non- degenerate for the functional equation ^ (and we denote 
{X, K) G ND{p)) if it satisfies the following conditions: 

• Spectral condition: the tangent space Tj^(^g)A4 has an invariant splitting for 
all 9 G T', 



(10) 



Tk{9)M - © © £K{e)^ 
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(12) 



(14) 



where S^f^gy ^Kie) '^'^'^ ^K{e) ^'"'^ stable, center and unstable invariant 
spaces respectively, i.e. 

DF(if(0))f;-« = 5--^). 

This splitting is analytic in 6. To this splitting we associate the projec- 
tions n^^g-j, 11/^(9) ^"■^ nK(e) fespectively, which are analytic with respect 

to e. 

Moreover, the splitting (|10|) is characterized by asymptotic growth con- 
ditions (co-cycles over T^): there exist < /ii,/i2 < 1 , M3 > 1 such that 
fiifis < 1, /i2M3 < 1 '^''^d Ch > such that for all n > 1 and 9 £ Dp 

\{DF) o K o T:'-1((?) X • • • X {DF) o K{9)v\ < Ch^^'M 
V e £K{e) 

\{DFy^ oKo Tj("-i)(6i) X • • ■ X {DF)'^ o K{6)v\ < Ch^i^\v\ 
<^=^ V E £k(0) 
and 

\{DF) oKo T^^^{0) X • • • X (DF) o K{d)v\ < Cft^a bl, 
(13) \{DF)'^ o K o Tj("-i)(0) X ... X (DF)-' o K{e)v\ < Chfi^\v\ 
<^=^ v G ^Kie)- 

• Furthermore, we assume that the dimension of the center subspace is 21. 

That is, the torus is as hyperbolic as allowed by the symplectic structure 
and there are no elliptic directions in the normal direction. 

• Twist condition: We introduce the notation 

N{e) = [DK{eyDK{0)]-\ 

p{e) = DK[e)N{9). 



Assume that the averages on T' of the matrices 

(15) Q{e)^DK{e + uj)'^J{K{e^oj))G{9) 
and 

(16) A{6) = P{6 + LuY \[DF{K)J{K)-^ P]{e) - [J{K)-^P]{e + uj) 
are non-singular. 

Remark 3.5. With a view to applications, we note that in Proposition l5.2[ we will 
show that we can deduce the existence of an invariant splitting from the existence of 
an approximately invariant one which satisfies the hyperbolicity Conditions (|lip - 
P^ . Consequently, Definition 13 .41 can be verified with a finite precision calculation 
on a given numerical approximation. We anticipate that the basic idea is that, if 
we can verify that for some operator B we have ||-B^|| < /i^ < 1 for some > 0, it 
follows that < C/i" for all n > 0. This gives a way to obtain all inequalities 

from finite computations. 

Remark 3.6. Note that since K is an embedding — hence DK{6) is one to one 
for aU 6* — and d > ; we have that DK{9)^DK{e) is invertible for all 9. 
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Remark 3.7. If wc take G{0) = J{Ko{0 + uj))-^DKo{9 + uj),Q becomes DK{e + 
ujyj{Ko{e + u))J{K{e + uj))-^DKQ{e + tj) « iV(e' + lo)-^ and hence one of the 
twist conditions becomes automatic because, under the smallness assumptions, 
avg (Q) := j^iQ{9)d9 k, avg(7V~^) and avg(7V~^) is invertible. Indeed, assume 
that V e Ker (avg {N^^)). Then u^avg {N~^)v — 0. This last expression is approx- 
imately 



which implies DK{9)v = for all G T'. Since DK{9) is one to one for all 9 
we obtain u = 0. Hence the condition on the invertibility of avg (Q) is just a 
quantitative statement of the fact that K is indeed an embedding. The condition 
on A is a twist condition. 

Remark 3.8. Note that if the torus K was exactly invariant (i.e. F oK = K oT^) 
then 



so that Conditions (fTT|) - ((T3)) are the usual growth conditions in the theory of nor- 
mally hyperbolic manifolds (see |Fen74| |HPS77|, lPes04| ). Of course, for our appli- 
cations, we only assume that the tori are approximately invariant. 

When the manifolds are Euclidean, the conditions (|lip - p^ make perfect sense. 
Nevertheless, if the phase space is a general manifold M., we have DF{K{9)) : 
Tk(0)M ^ TfoK{0)M. If FoK{9) ^ K{9+u}), then, we should write the conditions 
([TI])-([l3l) using connectors (see |HPPS70| 'l. 

We recall that 

Definition 3.9. A connector 5*1 is an isomorphism from TyA4 to T^M., defined 
when d(x,y) is small enough, such that — Id and S^Sy — S^, when both make 
sense. 

A concrete way of implementing the connectors is to take parallel transport along 
the shortest geodesic joining x, y (equivalently, the differential of the exponential 
map). 

In the case that we formulate the result in a general manifold, pT|) should be 
written 



(17) 

\{DF) o K o T:'\9)s2'^f_Q'^ X ... X (i?F) o K{9)v\ < C,^,^^\v\ 



and analogously the others. 

Remark 3.10. The technical reason why we introduced the extra parameter A in 
([7]) is the following: in the iteration of the KAM scheme, one has to prove that some 
equations are approximately solved up to a quadratic error. To this end, we have to 
show that some averages are quadratic in the error. To avoid these technicalities, 
we introduce this parameter A which allows us to cancel some terms in the equation 
so that we can reach the suitable approximate solution (See Propositions 14.181 and 
I4.19P . Then we use the exact symplecticness of the map to keep the parameter A 
under control. 




DFoK 



o T^-^{9) x---xDFo K{9) = DF" o K{9), 
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We can now state our main theorem, which provides the existence of a solution 
K to the functional equation ([6]) with F exact symplectic, provided we are given a 
sufficiently approximate one. 

Theorem 3.11. Let oj e D{k, v) for some n > >l. Assume that 

(1) F :IA G M. M is an exact symplectic map and U is an open connected 
set, which we will assume without loss of generality has a smooth boundary. 

(2) Kq G ND{po) (the embedding Kq is non- degenerate) in the sense that it 
satisfies the spectral condition in Definition \3.4\ and the average on T' of 
the matrices Qoid) and Ao{0) are non-singular, where Qq and Aq are as Q 
and A in Definition \3.4\ with K = Kq . 

(3) The map F is real analytic and it can be extended holomorphically to some 
complex neighborhood of the image under Kq of Dp^ : 

Br = {ze C^'^l 39 e {|Im6l| < po} s.t. \z - Ko{9)\ < r} , 

for some r > and such that |_F|p2(5^') is finite. 
Denote Eq — F o Kq — Kq o the initial error. Then there exists a constant 
C > depending on I, v, \F\c2(b^), \\DKq\\p,,, \\Nq\\p„, \\Ao\\p,,, |(avg (Ao))-^!, 
|(avg (Qo))~"'^|) |>^|ci(Sr) ihe norms of the projections ||n/('^'(0) llpo such that, if 
Eq satisfies the estimates 

(18) CkU-'"'\\Eo\\p, < 1 
and 

Cn^S-^-'WEoWp^Kr, 

where < S < min(l,po/12) is fixed, then there exists an embedding Koo S 
ND{pryo pq — 6(5) such that 

F o Koo = ° T^- 
Furthermore, we have the following estimate 

(19) lli^oo - KoWp^ < Ck^S-^'WEoWp,. 

Remark 3.12. The previous theorem provides a construction of whiskered tori 
without assuming the existence of action-angle variables for the original system. 
Moreover, the method of proof does not involve the sequence of transformations by 
symplectomorphisms, which is often used to prove this kind of results, but hard to 
implement numerically. 

Remark 3.13. It is important to remark that the non-degeneracy conditions we use 
in Theorem 13 . 1 II depend only on the approximate solution under consideration. As 
one can see, Definition l3.4l onlv depends on averages of the approximately computed 
solutions. This latter fact is useful in the validation of numerical computations. 
Indeed, numerical computations provide an approximate solution and this is the 
only information that is available. The non-degeneracy conditions needed to apply 
Theorem 13.111 can be verified by straightforward computations on the numerical 
approximation. 

This leads directly to the so-called small twist theorems. See Section [7.31 and in 
particular Proposition 17.11 and the subsequent comments for more details on the 
dependence of the constants on the non-degeneracy assumptions. 
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After introducing an additional term in the functional equation ([6|), namely 

FoK+ {J{Kn)-^DKn) oT^\ = KoT^ 

and performing a KAM iteration on [K, A) , the final task consists of proving that 
Aoo = using the geometry. This is done by using the exact symplecticness of F 
and a suitable representation of the center subspace. Indeed, the center subspace 
in T^(g)A^, which will be shown to be non-trivial, will be very close to the vector 
space spanned by DK{9) and its symplectically conjugate J [K {9))^^ D K {9) . 

3.3. Uniqueness. A natural question to ask is whether the embedding K provided 
by Theorem 13.111 is unique. Notice that if K is a solution of fB]), for any cr G R', 
KoTa is also a solution, hence one can only hope for uniqueness up to a composition 
with a translation on the right. 

The following theorem provides a local uniqueness result. We will see in the next 
section that there is a simple general argument that shows that uniqueness results 
allow us to deduce results for flows from results for diffeomorphisms. 

Theorem 3.14. Let F be exact symplectic and analytic in Br C A4. Let uj G 

D{k^v) for some n > 0,v > I. Assume Ki,K2 G ND{p) with p > are two 
solutions of equation ([6]) such that Ki{Dp) C B^, K2{Dp) C Br- Then there exists 
a constant C > depending on I, v, |F|c2(s^,), ||£>iiri||p, ||-/Vi||p, \ J\c^(Br)j II^iIIp? 
Iln^'^^g" lip, |(avg such that if for some t £ M' the norm \\Ki o — K2\\p 

satisfies 

(20) CK^p-^''\\KioTr-K2\\p<l 

with 6 = p/A, there exists a phase f G such that Ki oTf — K2 in Dp. Moreover 

\f-T\<CK^p-^''\\Kl-K2\\p. 

The proof of this theorem is postponed to Section [G] 

3.4. Result for flows. As a by-product of the previous uniqueness theorem, we 
get a result on the existence of invariant whiskered tori for flows. This follows from 
a time-one map argument (see |Dou82j ). The argument we present here comes from 
idlLGJV05[rCFdlL03j . 

Theorem 3.15. Let uj G D(k.,v) for some k > 0,i^ > I. Let {St)teM be the flow 
generated by a finite- dimensional analytic exact symplectic vector-field 

du 

Tt = 

where u : I ClR A4. Assume that there exists a time t = I and an embedding 
K G ND{p) for some p> such that Si o K{9) = K{9 + uj) for all 9 e TK Then 
for all time t G M, we have 

StoK{9) =K{9 + ujt). 

Proof. If we have Si o K{9) = K{9 + oj), then for all t this yields 

SioSto K{9) ^Sto {Si o K){9) ^Sto K{9 + uj). 

By Theorem 13.141 if || S't o X — iiTjlp is sufficiently small, which is achieved if t is 
sufficiently small, this implies that there exists a phase (f>{t) such that St o K{9) — 
K{9 + (f>{t)). From the flow property St+s = St o Sg and the fact that K is one to 
one, we have ^(i + s) — 0(i) -|-0(s). We now prove that the function (j) is continuous. 
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The map K from T' into its image is one to one and continuous over a compact (for 
the topology of T'). Then its inverse is continuous. This leads to the continuity of 
the function cf). Using this fact and the additivity condition we deduce, that for t 
small enough, (f){t) = (3t for some /3 e M'. Then in this case we have 

(21) StoK{e)^K{9 + l3t). 

Since both sides of ((2T|) are analytic with respect to t £ [0, 1] we obtain the result 
for all t e [0, 1]. Putting t = 1 we get (3 — uo. Expression ([2T|) shows that the torus 
if(T') is invariant by the flow. Since the torus is compact, the flow on it is defined 
for alH e M and hence holds for all t £ K. This ends the proof. □ 

In Section [9l we will give a more precise version of this result and a direct proof 
(i.e. a proof which does not pass through a reduction to a time-1 map). This 
is useful since the method of proof leads to numerical algorithms for differential 
equations. The direct proof can also be used as a model for results for some ill- 
posed partial differential equations (see [dlLS07| ). 

4. The linearized operator Dx^KJ^ui^, K) 

In this section, we describe the inductive step of the procedure. As most of the 
KAM proofs, it will be a modification of the classical Newton method. 
Using the Taylor theorem, given an approximate solution, we write 

+ A,K + A)^ K) + Dx,kT^{\, K){A, A) + 0(|(A, A)!^) 

and, following the idea of Newton's method, we look for (A, A) such that J-lj{X + 
A,K + A) is quadratically small so we are lead to consider the following equation 

(22) Dx,KrUX,K){A,A) = -E, 

where (A, K) is a pair satisfying approximately equation ([7]) with an error E{0) = 
J-i^{X, K){9) with 6 £ T'. Using the definition of the operator in ([5]), we see 
that the derivative of the operator can be written more explicitly as: 

DxmJ'U^ k){a, A){e) = G{e)A + DF{K{e))A{e) - A{e + lo). 

The study of the Newton equation is mainly done in three steps: 

• One projects equation ([22|l on the hyperbolic space and the center space, 
by using the invariant splitting (see Definition 13. 4p . 

• One reduces the equation of the projection on the center subspace to two 
classical small divisors equations. Thanks to a suitable change of coor- 
dinates on the tangent space (which does not use action-angle variables) 
these equations are then solved approximately (i.e. up to quadratic error) 
by using the extra variable A £ K' . 

• One solves (with "tame" estimates) the equations corresponding to the 
projections onto the stable and unstable invariant subspaces, by using the 
conditions on the co-cycles over T^^. 

Remark 4.1. We note that the equation on the center subspace will not be solved 
exactly. We will just solve it up to quadratic errors. The reason is that the change 
of variables mentioned in the above discussion will be constructed taking advantage 
of approximate identities obtained by differentiating with respect to 9 the equation 
for the initial error and applying geometric identities. The procedure of comparing 
the linearized Newton equation with the equations that appear taking derivatives 
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is very common in KAM theory. It is certainly used systematically in |Mos66bl 
IMos66a| IZeh76| . See [Zeh75l Section 5] for some remarks on the relation of these 
identities with a group structure of conjugacy problems. We note that some of 
these remarks in the above references work also for some semi-conjugacy problems. 

Of course, the above-mentioned strategy uses the non-degeneracy assumptions. 
In subsequent sections, we will show that these assumptions are changed only by a 
small amount, so that the procedure can be iterated. 

The main goal of this section is to prove the following result. 

Lemma 4.2. Consider the linearized equation 

(23) Dx^kTU^,K){A,A) = -E. 

Then there exists a constant C that depends on v, I, \\DK\\p, \\N\\p, ||n^'^'J'j ||p, 
IIGjlp, |(avg |(avg (Q))^"'^ I and the hyperbolicity constants such that assum- 

ing that S € (0,p/2) satisfies 

(24) Ck6-^''+'\\\E\1+\\G\\p\X\)<1 
we have 

(1) There exists an approximate solution (A, A) of (|23p . in the following sense: 
there exits a function E{9) such that (A, A) solves exactly 

Dx^KTU\K)iA,A) = -E + E, 

with the following estimates: for all S G (0,p/2) 

(25) l|A||p-25 < CK^r^-pii^, 

(26) \\DA\\p^2S<Ck'S-^--'\\E\\p, 

(27) \A\<C\\E\\p, 

(28) \\E\\p^s < Ck'S-^^''+''>\\E\\p\\TU\K)\\p. 

(2) // Ai and A2 solve the linearized equation in the previous approximate 
sense, then there exists a G M' such that for all S £ (0, p) 

(29) II Ai - A2 - DKi0)a\\p^s < C^^^C^^+D ||£;||p||.F„(A, K)\\p. 

Remark 4.3. The form of the previous inductive lemma corresponds very closely 
to Zehnder's implicit function theorem in |Zeh75| . Once Lemma 14.21 is proved, we 
then follow the strategy in [Zeh75| . The most crucial step is the verification of how 
the hypothesis of hyperbolicity are changed when the embedding changes in the 
iterative step. 

More precise information on the dependence of the constants C on the non- 
degeneracy conditions will be provided in Proposition 17.11 We anticipate that, 
roughly speaking , the constants C can be bounded by universal powers of the non- 
degeneracy constants. We postpone the precise formulation since it will involve 
some notations that will be developed along the proof. This power dependence on 
the constants has some applications to the study of tori close to resonance and to 
small twist theorems. 

We will need the following classical proposition (see |Riis76a| , [Riis76bj , [Riis75| , 
[dlLOlcj ) which provides existence of a solution together with estimates for small 
divisors equations. 
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Proposition 4.4. Let uj ^ D{k, v) and assume the mapping h ^ M. is analytic 
on Dp and has zero average. Then for any < a < p the difference equation 

v{9 + Lo) -v{d) = h{e) 

has a unique zero average solution v : ^ M., real analytic on Dp^„ for any 
< (7 < p. Moreover, we have the estimate 

(30) \\v\\p^^<CKa~''\\h\\p, 

where C only depends on v and the dimension of the torus I . 

Remark 4.5. It is important for our purposes to have estimates independent of 
the dimension of the manifold M since in a foUowup paper [FdlLSOS] we apply the 
procedure of this paper in an infinite-dimensional context. 

The independence of the estimates on the number of dimensions comes from the 
fact that we consider the supremum norm and the equation is solved component- 
wise. 

4.1. Geometric considerations. 

4.1.1. Isotropic character of the torus. We start by recalling the definition of isotropy. 

Definition 4.6. Let (A1,57) he a symplectic manifold. A submanifold M of Ai is 
isotropic if M C , where is the orthogonal space of Af with respect to the 
2- form D,. 

We formulate in our framework the well-known fact that a torus supporting an 
irrational rotation is isotropic. The manifold K{T^) is isotropic if the pull-back 
K*V,{9) vanishes for all 6* G T'. In other words, noting 

K*n{0){^,fj) ^<^,L{e)r^> 

for all ^, 77 G M', the isotropic character is equivalent to 

L{9) = DK{ey J{K{e))DK{9) = 

for all ^^ G T'. We first deal with the case of an exact solution of ([6]) (see Lemma 
14. 7p . The approximate case is the purpose of Lemma 14. 81 We note that the fact that 
exactly invariant tori are isotropic manifolds remains true for all irrational rotations 
and is well known |Zeh76j . The fact that approximately invariant tori carrying an 
irrational rotation are approximately isotropic seems to require that the rotation is 
Diophantine, see |dlLG JV05] . For the sake of completeness, we present the simple 
proofs of both results. 

Lemma 4.7. Assume that Ai is exact symplectic, K satisfies ^ and ui is rationally 
independent. Then L{9) is identically zero. 

Proof. Since F is symplectic we have 

F*D. = Q.. 

Consequently, this yields 

K*vt = K*F*vt ^{Ko T^yn. 

Since lo is rationally independent, T^^ is ergodic and this implies that K*n is con- 
stant and so is L{9). Using that Ai is exact symplectic, we have that K*il — dK*a 
and, the only constant form which is exact is zero. 

Similarly, a computation shows that L{d) has the form DLi{9)^ — DLi{9) for 
some matrix Li{9). Since the average on T' of DLi{9) is zero, we get the result. □ 
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Lemma 4.8. Assume that A4 is an exact symplectic manifold, F : Br ^ Ai is 
analytic and symplectic. Let K be real analytic on the complex strip Dp for some 
p > and such that K{Dp) C B^- Assume also that lu G D{k, v) and denote 

E = FoK + GX-KoT^. 

Then there exists a constant C depending on I, v, \\DK\\p, \F\(ji(^g^^, \J\c^(Br-) 
such that for all 6 G (0,p/2) we have 

(31) \\L\\p^,s<Ck5-(-^+'H\\E\\p+\\G\\p\M). 

Proof. We want to estimate the norm of the matrix L. RecalUng F*il — il, one 
gets 

K*fi ~{Ko T^yn = E*n - {Gxyn. 

Performing the same computations as in |dlLG JV05] . this leads to the following 
equation 

L - LoT^ = g, 

where g is a function on T' such that (here we just use Cauchy estimates) 
\\g\\p-5<C5-\\\E\\p + \\GU\\). 
We now make use of Proposition 14.41 to complete the proof. □ 

Recall that we are assuming that K is an embedding. Hence the range of DK 
is ^-dimensional. 

4.1.2. Vanishing lemma. This section is devoted to an estimate which allows to 
control the extra parameter A through the iterative step. We consider the functional 
equation 

FoK + GX = KoT^ + E, 

where F is exact symplectic (see Definition I3.2[) and G = [J{Kq)~'^DKq\ oT^. 

Recall that A € and Kq G ND{po). Note that the term {J{Ko)~^DKo) o 
is very close to { J{K)^^ DK)oT^ and hence close to the center subspace associated 
to the torus KiJ^). 

The following lemma provides the desired vanishing result. 

Lemma 4.9. Assume F maps A4 into itself and lo G D{k,, v). Let K G ND{p) be 
a solution of 

(32) F o K + GX^ K oT^ + E, 
with G — [J{Kq)^^ DKq] o and X is such that 

ll^llp+l|G||p|A| <r, 
^ ' \\K-Ko\\p<r, \\DK -DKo\\p<r, 

where r > is sufficiently small (precise conditions will be given along the proof). 
Assume furthermore that 

(1) F is exact symplectic. 

(2) F extends analytically to a neighborhood of K{T^). 
Then, there exists a constant C such that 

\M<C\\E\\p. 
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Proof. Wc follow a method used in [JdlLZ99] . We refer the reader to the Figure [T] 
for an illustration of the method. 
We denote by 

(34) = eT'-i 

and similarly uji = {oji, . . . ,uji^i,uji+i, . . . ,uji) e MJ'^-^. 
We also denote a-g, : T ^ T' the path given by 

(35) a^,-^(7/) = (0i,...,0,_i,ry,0,+i,...,0O- 

We will compute Jji-i J^oa. F*ct in two different ways. On one hand, using the 
fact that F is exact symplectic, we have 

(36) / F*a^ I a + dW^ 

JKoa, 

On the other hand, using ([3^ 



(37) / F*a = / a = / a + R,, 

JKoa.j, JFoKoa.g, J [KoT^-G\)Q(y. g , 

where \R,\ < C\\E\\p. 

Since we want to compare the last integrals in p6p and ([37]) it is natural to intro- 
duce a two-cell whose boundary is the difference between the two paths K oa,- g,^^_ 
and {K oT^~ GX) o • We denote B^^ this two-cell, which we parametrize by 
(C, ri) e (0, 1) X (0, 1) as' follows: 

By Stokes's theorem, since da — il, we have 



(38) / a = / a+ VL. 

We have 

/ n= f /■'l7s,,.,(c,,)(9ci?,,,-,.,(C,r?),9,i3,«.,,(e,r,))dCrfry. 
Note that 

and 

Using the previous expressions 
"s.,^^,,(C,r,)(9cS^,e.,A(^'^/)'^';^«J.,A(C'^)) 

= - A^Go fT,^^^+-^(r7)^J(i?,,,-^,,(e,r;))(9,.X o (t?) - defi o (7?)A0. 

Using ([33)) . we have 

•^(^.,e..A(^' '?)) - o <T,^o.^^^{v)) + 0(|A|) 

= J(i^ooa^,-^^.^(77))-hO(r)+0(|A|). 
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Therefore, we end up with (using the expression of G{9) — J{Kq{0)) ^DKq{9)) 
[ n = - f f \'^[DK^J{Ko)-^J{Ko)de^Ko]oa^^,^.Xr^)d^dr, 

JB.g,^ Jo Jo 

+ 0{r\X\)+Oi\\n + Om\p)- 
Joining these expressions for all values of i and integrating over T'^^ we get 



(39) / / n = x 



Q 



+ 0(r|A|)+0(|An + 0(||£;||,) 



where Q = DKj DKq. Since DKq has rank / then the matrix Q has rank I. See 
Remark 13.61 

We now integrate with respect to di both ([36|) and ((37|) . By a simple change of 
variables we have that the following integrals are equal 

d9^ / F*a = / de, F*a. 

Therefore, from (|38p we obtain 



Now, equation (f39|) . the fact that Q is invertible, the assumption 

||£;||p + ||G||p|A| <C, 

and r sufficiently small (this is the condition we imposed in (|33p ). leads to the 
desired result invoking the implicit function theorem. □ 

Remark 4.10. The assumption in Lemma l479l that 



||£;||p + ||G||p|A| <C 

will be an inductive assumption in the iteration of the KAM method that we will 
deal with later. 

Remark 4.11. In the KAM iteration, we will generate a sequence {A^, -^n}^^]^ of 
approximations of the solution (Aoo, ^oo) of the equation 

F o K + G\^ K oT^. 

As a corollary of Lemma l479l the sequence {A„}„gpj in the KAM iteration converges 
to since ||-En||p„ converges to 0. 



4.1.3. Basis of £'jf{e) "when K is an exact parameterization. To avoid the use of 
action-angle variables we are going to perform a change of variables, using the 
geometric structure on the tangent bundle. 

For that we will first find a useful basis of the center space f ^(g-j in the case that 
K -.T' ^ Misa. solution oi F o K = K oT^. 

Here we are assuming that the dimension of the center subspace is 21 and hence 
^K(e) [dlLGJVOS] . the authors studied the case when d = I, i.e. the 

dimension of the range of K{T'') is half the dimension of the space and the tori are 
Lagrangian submanifolds. 
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Figure 1 . Illustration of the vanishing lemma 



In |dlLGJV05] it is shown that, in the Lagrangian case, the perturbative equa- 
tions can be studied very conveniently applying the change of variables given by 
the following matrix 

(40) [DK{9),J{K{e))-'^DK{9)N{e)]. 

In the Lagrangian case, the range of ([^0]) is the tangent space of the manifold at 
K{9). In our case, however, the range of (|40)) is not the whole space, but it will be 
a very good approximation of the center space. Then, we can apply a method very 
similar to the method in [dlLGJVOS] for the equations in the center directions. The 
hyperbolic directions will be solved by other methods. 

By the symplecticness of F and the dynamical properties we have that the matrix 
of the symplectic structure with respect to the splitting f^^g-, © ^xie) ® ^k(9) 
the form 










JSU 


(41) 


j{Kie)) = 


JCC 







\jus 









where J^'^ is an antisymmetric form and J'^^{es,eu) — — J"^(e„,es). 
Indeed from 

u^J{K{9))v = nK^g){u,v) = np^^K(9)){DF''{K{9))u,DF^{K(9))v), n £ Z, 

we deduce, sending n +cxd and using the hyperbolic conditions (expansion/contraction 
properties), that J{K{9))v = in the following cases 

• w G U and V £ 

• we and V e 5^(3) U El^g^ 
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which imphes the form (j4T]). See also [DdlLS08| . The form jH]) proves that 
J{K{9))~^ sends the center subspace into itself. 

Since range!) if (6*) is the tangent space of the torus if(T') and the dynamics on 
the torus is conjugated to a rotation, DK{9yM} is contained in S'^^gy Moreover the 
previous property of J{K)^^ implies that J{K{9))^^DK{6)M} also is contained in 
£:^(g). Instead of J{K{9))-^DK{0) we will consider the matrix J{K{0))-^DK{e)N{0) 
where N{9) is the normalization / x Z-matrix N{9) = [DK{9)^ DK{9)]~^ introduced 
in (jl4p . Both have the same range because N{9) is non-singular. The role of N is 
to provide some normalization for the symplectic conjugate. 

Now we check that the range of [DK{9), J{K{9))-'^ DK{9)N{9)] is 2/-dimensional. 
Indeed, assume that there is a linear combination 
I I 
f = Y,a,DK{9)e, f^jJiKi9))-^ DK{9)N{9)ej = 0. 

Then, for 1 < fc < /, using the isotropic character of Tx(9)K{T'') 

I 

= n{DK{9)ek, f) = J2 Pi4DK{9y J{K{9))J{K{9))-^DK{9)N{9)e, 
i=i 

This calculation shows that / reduces to X]j=i '^jDK{9)ej. Moreover, for 1 < fc < Z 

- n{J{K{9))-^DK{9)N{9)ekJ) 
I 

= J2 a,elN{9yDK{9yj{K{9))-^J{K{9))DK{9)e, 
I 

Hence aj = f3j = for all j = We conclude that 

i-ange[DK{9), J(K{9))-^ DK{9)N{9)] = E'^^^y 
Finally we define 

(42) M{9) = [DK{9),J{K{9))-'^DK{9)N{9)]. 

4.2. Solving the linearized equation on the center subspace. This section is 
devoted to the study of equation projected on the center subspace. We denote 

Projecting the linearized equation (|23p into the center space, we end up with the 
following equation 

(43) n^(,+^)G((?)A + DF[K{9))^'^{9) - A=(0 + u:) = -£^(9), 

where E^{9) = Il'^K(e+^^m 

In Section 14.2.11 we will develop several identities and approximate identities 
that have a geometric nature. These identities will be used to reduce the equation 
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on the center to constant coefficients equations of tlie form considered in Propo- 
sition |431 One important step is accomplished in Section [4.2.21 where we use the 
geometric identities and the theory of hyperbolic systems to obtain an approximate 
representation of the center space. Once this material is developed, we can establish 
the main result of this Section, Proposition 14. 191 

4.2.1. Normalization procedure. In the following, we construct a suitable represen- 
tation for the matrix DF{K{9))M{d). Recall that the 2d x 2/-matrix M is given 

by 

(44) M = [DK, J{K)-^DKN]. 

As a motivation, we first consider the case when K \s a. solution of ([5]). We 
search for a matrix S{0) satisfying 

(45) DF{K{e))M{0) = AI{0 + uj)S{d), 

where S{9) is upper triangular with identity matrices on the diagonal. Explicit 
expressions for S will be given later. 

Differentiating equation © with respect to 61, we get 

DF{K{e))DK{e) = DK{e + uj). 

This shows that S{6) has the form 

'Id, A{ey 
,0/ B{e)^ 

where A{6) and B{d) are / x I matrices. We will see that the choice of the second 
column of M and the symplectic structure forces that B{9) — Id;. Then, it will be 
easy to compute an expression for A. 

Indeed, from (l44|), jUl) and ([46]) wc should have 

(47) [DF{K)J{K)-'^DKN]{e) = DK{e + uj)A{e) + [JiKy^ DK N]{e + uj)B{e). 
By the isotropic character of X(T') we have DK^ J{K)DK = 0. Hence 

(48) [DK^ J{K)]{9 + Lj)[DF{K)J{K)-^DK N]{e) = [DK^ DK N]{9 + Lj)B{e). 
Also by the symplecticness of F 

J{K{e + uj))DF{K{9)) = J{F{K{e)))DF{K{e)) = [DF{K)-^ J{K)]{e). 
Then the left-hand side of (H51) becomes 



(46) 



DK^ {e + uj)[DF{K)-^ DK N]{9) = [DK^ DK N]{e) = Id;. 

With this we conclude that B{9) = Id;. 

To obtain the expression of A{9) we multiply (gll) by {DK N){9 + w)^. Using 
=N and NDK'^DK = Id; we get 

(49) A{e) = P{9 + \[DF{K)J{K)-'^P]{e) - [J{K)-^P]{e + uj) 

We sum up the previous computations in the following lemma. 

Lemma 4.12. Let K he a solution of equation Then we can write 

DF{K{e))M{9) = M{e + uj)S{e), 

with 
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A{e) - p{e + \^DF{K)j{K)-^p]{e) - [j{K)-^p]{e + w) 

where the notation P{9) = DK{d)N{9) was introduced in (|14p . 

The matrix M(9) is not invertible since it is not square. However we can derive 
a generalized inverse for M{9). As a motivation for subsequent developments, we 
first present Lemma [4.141 which deals with the geometric cancellations in the case 
of an exactly invariant torus. The case of interest for a KAM algorithm — when 
the torus is only approximately invariant — will be studied in Lemma 14.151 as a 
perturbation of Lemma 14.141 

A straighforward calculation shows that 

L Idi 

-Idi DK^ J{K)-^ DK N J 



(51) M^J{K)M = 

Remark 4.13. When = — Id2d we have J^^ = — = J and then 



(^2) ^^"^W^^=(-H, N-LN 

If moreover if is a solution oi F o K = K oT^ the right-hand side matrix of ([52 

Id, 



reduces to the standard symplectic matrix Jo 



Id, 



Lemma 4.14. Let K he a solution of ([6]). Then the matrix J{K)M is invertible 
and 



N^DK^J{K)-^DKN -Id, 
Id, 

Proof. It follows from (jSTjl and the isotropic character of the invariant torus, i.e. 
L = 0. □ 

We now establish a similar result for approximate solutions, i.e. solutions of ([7]) 
up to error E{0) = Tuj{X, K){9). We can expect this type of normalization to be 
true if the error and A are small enough. Following the calculations in Lemma |4T2l 
we obtain: 



0{E'',DE''). 



More precisely, we introduce 

(53) e{e) = DF(K{0))M{0) - M{e + oj)S{e), 

where S is given by ([50)1 . If wc denote e[9) — (ei(^?), 62(6*)), a simple algebraic 
computation yields 

ei{e) = DE^e) - DgG''{e)\ 

62(6*) = [{DFJ-^){K)DK N]{e) - DK{e + uj)A{e) - [J-^ DK N]{e + uj) = 0{E,DE) 

by the choice of A, where G^(6i) = n^(e+^)G(6i). 

The next step is to ensure the invertibility of the 21 x 2Z-matrix J{K)M . 
According to expression ([5T|) . we can write 

M{eY j{K{0))M{e) = v{e) + R{e), 
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where 

Idi 



V 
and 

R 



-Ml N ' DK ' JiK)- ' DK N 



L 0^ 

Oy 

We have the foUowing lemma, providing the desired invertibility result under a 
smallness assumption on E, namely ((54|) in the next lemma. Note that (|54p has 
the same form as (f24| , but the constants could be slighly different since ([M)) should 
also accomodate (I33p . which is implied by conditions of the same form. 

Lemma 4.15. There exists a constant C > such that if 

(54) CkS-^''+'^^\\E\\p< 1/2 

for some < S < p/2 then the matrix {9)J{K{9))M{9) is invertible for 9 G 
Dp-2S and there exists a matrix V{9) such that 

{M{fff J{K{9))M{9))-^ = V{9)-^ + V{9) 

with 

V{9)=(j2^V{9)-^R{9)t)v{9)-\ 

k=l 

where the series is absolutely convergent. Furthermore, we have the estimate 

(55) \\V\\p-25<C'kS-^^+'^\\E\\p, 

where the constant C > depends on I, v, |-F|ci(Sr)i l"^lci(Sr); ll^-^^llp; ll^llp and 
l|n^(e)llp- 

Proof. The matrix V{9) is invertible with 
V 

We can write 



1 _ (N^DK^J{K)-^DKN -Id, 
Id, 



M{9)'^J{K{9))M{9) = V{e){ld2i + V{9)-^R{9)). 

To apply the Neumann series (and consequently justify the existence of the inverse 
of Id2, + V^'^R as well as the estimates for its size), we have to estimate the term 
V~^R. According to Lemma we have the estimate for L 

\\L\\p^28<Cn5-^'^+'H\\E\\p+\\GU\\) 

for all 5 e (0,p/2). Using Lemma l4!9l this leads to the estimate 

\\V-^R\\p-25<Ck5-^-+^^\\E\\p 

for < (5 < p/2, where C > depends on v, \ J\c^{b^), \\DK\\p, \\N\\p 

and ||n^^g^||p. Because of assumption (|54)) . we have that the right-hand side of the 
last equation is less than 1/2. 

Then the matrix Id2, + V [9)^^ R{9) is invertible with 

This ends the proof of Lemma 14.151 □ 
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4.2.2. Identification of the center space. In this section, we identify the center space 
as being very close (up to terms that can be bounded by the error) to the range of 
the matrix M introduced in ([l^ . see Proposition 14.161 This will allow us to use 
the range of M in place of S'^^g-^ without changing the quadratic character of the 
method. 

Proposition 4.16. Denote by I^Kie) the range of M{9) and by 11^,^^^^ the projection 
onto ^Kie) according to the splitting Sx(e) ® K(e) ® ^K{e) ■ 
Then there exists a constant C > such that if 

5~'\\E\\,<C 

we have the estimates (here distp stands for the distance between subspaces at the 
Grassmanian level) 

distp_25(r^(fl),£^(e)) < C6-'\\E\\, 

m'k(s)-nl^s)\\p-2s<cs-^E\i 

for every S G (0,p/2) and where C, as usual, depends on the non-degeneracy con- 
stants of the problem. 

Proof. Of course, the two inequalities in (|56p are equivalent. 
From (|53p and Cauchy estimates, we have: 



d:istp^s{{DFoK)TK(e),TK(e)oT^) < C5-'\\E 



IP- 



Using again equation (I53p and iterating it, we obtain for n > 1 
DF{K{e+nuj)) X • • • X DF{K{e))M{0) = M{0+nuj)S{e+{n-l)uj) x • • • x5(6l)+i?„, 
where 

\\Rn\\p-5<Cn5-^\\E\\p 

and Cn depends on n. 

Since S{9) is upper triangular with Id; on the diagonal, we have: 



S{e + {n-l)uj) x ••■ y.S{e) 



Id, A{e + {n - + ■ ■ ■ + A{e) 
Idz 



Therefore, by induction, we have for every rt G N 

\\DF[K{e + nuj)) ■ ■ ■ DF{K{e))M{e)\\p^s < Cn + Cn5~^\\E\\p. 

Identical calculations give that 

\\DF~\K{e - uLo)) ■ ■ ■ DF-^[K{e))M{e)\\p-s <Cn + C^5-^\\E\\p. 

Note that, given any /ia > 1 (as in Definition [33]), there exists an integer n^j > 
such that for all n > n^g, we have Cn < /ig . Consequently, choosing such n^j there 
exists a constant C such that if the error satisfies 

S-^\\E\\p<C 

we have Cn + C„(5~-^||i?||p < /Xg. In other words, the above estimates hold for 
all sufficiently large n, provided that we impose a suitable smallness condition on 
5-'\\E\\p. 

As a consequence, ^K{e) is an approximately invariant bundle, and we also have 
bounds on the rate of growth of the co-cycle both in positive and negative times. 
Using standard tools in the theory of hyperboHc systems (see Proposition IS . 21 below 
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where we prove the resuh for all the bundles), this shows that indeed one can find 
a true invariant subspace £K{e) close to r^(9). Since this invariant subspace should 
be of the same dimension of the center space ■S'^^g), we deduce that 

£k{9) = £k{9)- 

See also Remark 15.61 below. □ 



4.2.3. Final estimates of the solution on the center subspace. We can now finish the 
solution of equation ((23l) on the center subspace. We recall the linearized equation 
around (A, K) projected on the center subspace: 

(57) G"(6')A + DF{K{9))A%e) - A^id + u) = -E^iO). 
We make the change the unknowns in ((57|) 

(58) A^ie) = M{e)w{e) + e{e)w{e), 

where 

which was estimated in Proposition 14. 161 
Substituting ((58| into equation ((57|) we get 

(60) DF{K{9))M{e)W{9) - M{9 + uj)W{9 + uj) 

(61) = -£""{9) - G'{9)K + e{9 + u)W{9 + u) - DF{K{9))e{9)W{9). 
We anticipate that the term eW will be quadratic in the error. Similarly, writing 

we also anticipate that the term eGA will be quadratic in the error. Since the 
function G will be chosen to be J{Kq)~^DKq o T^, namely in Txg(^g^^^'j, we drop 
the index from G*^, writing G directly. As a consequence, we will ignore these two 
terms and consider instead the equation 

(62) DF{K{9))M{9)W{9) ~ M{9 + lu)W{9 + uj) = -E''{9) - G{9)A 

which differs from the linearized equation in the term (eW) oT^ — DF{K){eW) — 
eGA. Note that, ignoring this term we obtain an equation where all the terms are 
in the range of M . 

We muhiply equation ([62]) by M{9 + uj)^ J{K{9 + uo)) 

\M^J{K)]{9 + u)DF{K{9))M{9)W{9) - [M^ J{K)]{9 + u)M{9 + lu)W{9 + lu) 

^ - [M^ J{K)]{9 + uj)[E''{9) + G{9)K]. 
Using Lemma [4. 151 (invertibilitv of J{K)M) and equation (|53|). we can write 



(63) 



Id; A{9)\ 
Oi Id, ] + 



W{9) - W{9 + Lu)= pi{9) +P2{9) 
- [M^J{K)M] {9 + u)-^[M'^ J{K)] (9 + uj)G{9)A, 

where 

(64) B{9) = [M^J{K)M]{9 + uj)-'^[M^J{K)]{9 + uj)e{9), 

(65) pi{9) = -V{9 + uj)-\M^J{K)]{9 + u)E''{9) 
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and 

(66) P2i9) = -¥{9 + ij)[M^J{K)]{9 + lo)E^{9). 

In the following lemma, we sum up the previous computations and estimate the 
terms in equation 



Lemma 4.17. Assume u) G D{k,,v) and 5 and \\E\\p satisfy ([54|l . Equation ([62 
can he written in the form 



(67) 



Id, A{9) 
Oi Id, 



B{9) 



W{9) - W{e + Lu)=pi (9) + p2 (0) 



iM^ J{K)M]-\9 + uj)[M^ J{K)]ie + uj)G{9)A, 



where the matrix B and the vectors pi and P2 are given by expressions (j64p . (j65p 
and (j66p respectively. 

The following estimates hold: 

(68) lbillp<ci|£;i|p, 

where C only depends on |J|ci(_Br)' ll-^llp' ll-^-^'llp '^'^'^ l|n^(9)||p- For p2 and B 
we have 

(69) \\P2t-2,<CK5-^'^+^)\\E\\l 
and 

(70) \\B\l^25<C5-W\E\\, + \\\), 

where C depends I, v, \\N\\p, \\DK\\p, \J\c^{b^) and \\W^f^g-^\\p. 

Proof. Since the matrix V^^ does not depend on L the estimate (1551) is obvious 
from the formula (|65p for pi{9). 

According to the proof of Lemma [4.15l the estimate ([M|) then comes from estimate 
(l55t. We turn to the estimate on B. We have 



B{0) = {V{9 + uj)-'^ + V{9 + uj))M{9 + uj^ J{K{9 + u))e{9). 
This leads to 

\\B\\p-2S <\\V{9 + u)-^\\p-2s\\M{9 + uoY J{K{9 + u))e{9)\\p.25 
+ \\V{9 + Lu)M{9 + Lu)'^J{K{9 + ujj)e{9)\\p^2S- 
Therefore, using estimate (j55p and Cauchy estimates, we end up with 

\\B\\p-25 < C{d-'\\E\\p + 5-i|A| + KS-^-'+'^WEWpiS-'WEWp + S-'\X\)). 
This leads to the desired result thanks to the smallness assumption on ||i?||p. □ 

4.2.4. Approximate solvability of the equations on the center sub space. This section 
is devoted to solving approximately (up to quadratic error) the linearized equation 
(|67)) . as is usual in KAM theory. 

To this end, we introduce the following operator 

CW{9) = {^"^'^ W{9) - W{9 + .;). 



28 



E. FONTICH, R. DE LA LEAVE, AND Y. SIRE 



Equation ()67|) can be written as 

(71) cw{e) + B{e)w{e) =pi{e)+p2{e) 

~ [m'^ j{K)M]{e + uj)-^[m'^ j{K)]{e + uj)G{e)k. 

According to estimates in Lemma [4. 171 we have p2 — 0(||i?||p), pi — 0{\\E\\p) and 
B = 0(||_E||p + |A|). Solving approximately equation ([7T|) with an error "quadratic" 
in E does not affect the convergence of the Newton scheme. See jZeh75| for an 
abstract discussion and |Zeh76| for several concrete applications. 

Equation (|7l]) does not fit into the framework of Proposition H^ since the average 
of the right-hand side is generically non-zero. However, by using the increment 
parameter A, we can make this average equal to zero. Furthermore, equation (j7ip 
has two unknowns (the two symplectic coordinates). Thanks to Lemma |4.15[ one 
can write the term 

[M'^ J{K)M]-^{e + uj)[M'^ J{K)]{e + uo)G{e)A ^ qi{e)A + q2{0)A, 

where the matrix qi (which is 21 x /) is 

qi{e) = V{9 + u)-^M{e + ujyj{K{9 + u))G{e) 

and q2 satisfies for all 5 e (0,p/2) 

\\q2\\p-2S<CK5-^''+^'>\\G\\p\\E\\p, 

where the constant G depends on /, v, ||A^||p, H-DXHp, \F\c^{b,^)^ \J\c'^(Br) 
l|n^(e)llp- 

We define an approximate solution of (j7ip as a solution of the following equation 
(1721) . obtained by removing the terms containing B and q2 from the equation (|67p . 
which was equivalent to ([62|) . We recall that (j62|) was obtained from the Newton 
step by removing the terms that contained e. As we will see, all these eliminations 
do not change the quadratic convergence of the method. Consider now 

(72) Cv{e)=pi{e)^q^{e)k. 

Thanks to the non-degeneracy conditions (see Definition 13. 4|) . we obtain the 
following result. 

Proposition 4.18. Assume lo e I?(k, v) and (A, if) is a non- degenerate pair (i.e. 
(A, A') G ND{p)). If the error \\E\\p satisfies (|54p and the smallness assumptions 
in proposition \4.16\ there exist a mapping v, analytic on Dp_2S o.'^d a vector A G M' 
solving equation (j72p . 

Moreover there exists a constant C > depending on v, I, ||A'||p, |(avg (Q))^"'^], 
(avg(A))~-^|,, ||A^||p and ||n^jg^||p such that 

\\v\\p-28<CK^b-^''\\E\\p 

and 

\t.\<G\\E\\p. 

Proof. We denote R{0) the right-hand side of equation ([72|) . i.e. we solve 

(73) Cv{9) = R{e), 
with 

R^pi- qiA. 
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We now decompose equation ()73p into symplectically conjugate coordinates, i.e. 
V — {vi,V2)^ , R{0) = [Ri{9),R2{9)Y . Therefore, equation ([73]) is equivalent to 

vx{e) + A{e)v2{e) ^vx{e + Lo) + Ri{0), 

V2{e)=V2{9 + Uj)+R2{9). 

A simple computation shows that 

R2{9) = -[DK^J{K)]oT^G{e) 
We choose A G M' such that 

avg(i?2) =0. 

According to Proposition l4.41 if avg {R2) — the equation in V2 admits an analytic 
solution with arbitrary average on Dp^s and we have the estimate 

(74) \\v2\\p-5 < Ck6--'\\R2\\p + \SiYgiv2)\. 

Then we choose avg(u2) such that avg — AV2) = 0, which allows us to solve 
uniquely the equation in vi, the function vi being of zero average. Furthermore, 
we have the estimate 

\\vi\\p-2S < CkS-'WRi - Av2\\p-s. 

We now turn to the estimates. First we estimate A. The vector A G is such 
that 

avg (^DK^iu; + e)J{K{uj + e)){E%e) + G(6')A)) = 0. 

This leads to 

avg (^{DK^iuj + e)J{K{uj + 9))G{9)'^ A 

= - avg {{DK^{u + 9)J{K{u + 9))E''{9)^ . 

Note that by the definition of P and the fact that N is symmetric, the matrix which 
applies to A is the average of Q which, by hypothesis, is invertible. This leads to 
the desired estimate for A. 

We now estimate the solution v. From the expression of R and the value of A 
obtained above, we have that there exists a constant C such that 

\mp<c\\E\\p, 

for i = 1,2. Furthermore, we choose avg (1)2) such that avg (i?i — ^^2) = 0, i.e. 

avg('U2) = avg(A)~i(avg(i?i) - avg(Ai;^)), 

where V2 = V2 + avg (172). This is possible since by the twist condition avg (A) is 
invertible. Thanks to estimate (|7^ . this leads to the desired result. □ 

We now come back to the solutions of (|43)) . The above procedure allows us to 
prove the following proposition, providing an approximate solution of the projection 
of Dx,kTuj{X, K){A, A) = —E on the center subspace. 

Proposition 4.19. Let (A, W) be as in Propo.sition \4-. 1^ and assume the hypotheses 
of that proposition hold. Define A''{9) = M {9)W (9) + e{9)W (9) and obtain W and 
A as indicated above. 

Then, (A, A^) is an approximate solution of (|43|) and we have the following 
estimates 

\mp-2S<C^6-'-^\\E\\p, 
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|A|<C|li?||p, 

where the constant C depends on v, I, |(avg (Q))""'" |, |(avg (A))^^ |, ||-/V||p, 
and \\^x{e)\\p- Moreover 

(75) 

pA,K^c.(A,if)(A, A^) + E'\\p^2S < C«;3r(3.+i)(||^||2 ^ ii^ii^i^i) ^ cS-'+''\\E\\l 

where the constant C depends on I, k, v, |i^|ci(Sr)j II^-^IIp; [(^"^g 
|(avg(Q))-i| a«d||G||p. 

Proof. The first estimate comes from tfie previous Proposition 14.181 
For the second one (|75p. we recall that we have: 

Dx.,KT^{K,^''){e) + E%e) 

(76) ^-[M'^j{K)M]{e + uj)[B{e)v(e)-p2(e)]-q2{e)K 

+ e{e + uj)w{e + uj)- DF{K{e))e{e)w{e). 

The first term in the right-hand side of ([7^ is estimated in Proposition l4.18l see 
estimates ([69|) - ([70|) . The second one comes from the vanishing Lemma [4.91 The 
third term is estimated in Proposition 14.161 □ 

4.3. Solving the linearized equations on the hyperbolic subspaces. Ac- 
cording to the splitting p^ . there exist projections on the linear spaces and 
^^(9)- '^^^ analytic regularity of the splitting implies that the dependence of these 
projections in 9 is analytic in the same domain as the spaces. We denote n^^^^^^^ 
(resp. n^jg^^j) the projections (of base K{9 + lo)) on the stable (resp. unstable) 
invariant subspace. 

We project equation ((22|) on the stable and unstable spaces to obtain 

(77) n^(,+^) (g(0)A + DF{K{e))m - A(e + 0,)) = -n^K(e+u)E{e), 

(78) n^(,+^) (G{e)K + DF{K{e))m - A(0 + u)) - -nl^e+^)E{e). 

Furthermore, thanks to the invariance of the splitting, we can write 
^'Kie+^)DF{K{e))m = DF{K{e))Wj,^,)m 

for the stable part and 

^l(e+.)DF{K{e))m = DFiK{9))ni^,^Ai0) 

for the unstable one. Introducing the change of variables 9' — T^{9) and the 
notation A'*'"(6'') = n^"^,-^ A(6''), equations can be written in the following 

form 

(79) DF{K) o r_^(0')A^(T_^(0')) " A^(0') = ~E'{9',K), 
where 

^^(0',A) = n^(,,)(G(T_^(0'))A) +n^(,,)i?or_^(^?') 

and 

(80) DF{K) oT_^{9')/S}\T_^{9')) - A"(6i') = -£;"(6i',A), 
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where 

^"(0',A) =n^.(,,)(G(T_40'))A) +n^(,,)i?oT_40')- 
The following proposition provides an existence result together with estimates 
for equations ([75)) - ([50)1 . 

Proposition 4.20. Fix p > 0. Then equation (17^ (resp. ([5(1)1 ] admits a unique 
analytic solution A'' : Dp — > ^K(e) (''"^^P- ■ ~^ ^K(e))' Furthermore there 
exists a constant C such that 

(81) ||A^'"|lp<C(||i?||p + |A|), 

where the constant C depends on the hyperbolicity constant fii (resp. (12), the norm 
of the projection ||n^^g-j||p (resp. \\^fc{e)\\p) ll^(^)llp and the constant Ch involved 
in m\> (resp. fll)). 

Proof. We give the proof for the stable case, the unstable one being similar and left 
to the reader. Using equation ()79p iteratively, we claim that its solution is given by 

00 

(82) A^(^?')=E(^W°7^--(^?')x---xi?F(i^)or_^(0'))^'(r-(fc+i).(0'),A). 

k=0 

Using the condition on the co-cycles over T^^ (see equation ([TT))). the series con- 
verges uniformly on Dp and one can estimate 

oc 

(83) \\A^\\p<c,,\\E^\\pY.^,'l<cm\p + \M) 

k=0 

since fii < 1. Once we know that the series converges uniformly, we can rearrange 
the terms and get that ([82)1 is indeed a solution. The proof in the case of the unstable 
space follows in the same way, multiplying equation ([80)1 by {DF{K) oT^^)~^ and 
using the Condition ([T2)l on the co-cycles. □ 

5. Iteration of the Newton step and convergence 

In the following we describe precisely the iteration of the Newton method. As it 
is standard in KAM theory, we show that if the initial error ||-Eo||po small enough, 
one can choose the domain loss, so that the iterative scheme converges to a solution 
of ((7]l which moreover is close to the initial approximate solution. As a consequence 
of the vanishing lemma (i.e. Lemma one gets A = and then a solution of 

F oK = K oT^. 

In the rest of this section, we are under the assumptions of Theorem 13.111 

5.1. Estimates for one step of the Newton method. Recall that we have 
implemented a step showing that, given an approximate solution, (A„j_i, K^-i) of 
([7|) , which is non-degenerate in the sense of Definition l3.4l and satisfies the conditions 
([33)1 of Lemma 14.91 and ([24)1 of Lemma 14.21 , then we find an approximate solution 
(A,„^i, Am_i) of the Newton equation. That is, we can find 

with Em.^i{0) = .7\j(A„i-i, Km-i){d) and Rm "quadratically" small. If E is defined 
in Dp^ -^, the Newton correction Am_i is defined in a smaller domain Dp^, pm = 
Pm-i — Sm- The precise results on the step are collected in Lemma 14.21 and the 
description of the step is given along the proof. 
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The next result Proposition 15. H makes precise the observation that, if we can 
define F o Km then it is possible to show that the new remainder is quadratic. 
Furthermore, we will show that the change in the non-degeneracy assumptions can 
be estimated by the size of the error. 

The assumption that F o Km can be defined, requires only that the range of 
Km = Km-i + Am-i docs not get out the domain of F. This will be implied 
by smallness assumptions on A that, using the conclusions of Proposition 14.21 are 
implied by assumption ([55]) . As it will turn out, the assumption (|85p is stronger 
than (211) so that ([55)) is enough to ensure that we can carry out a Newton step as 
indicated. 

In subsequent sections, we will show that if we choose the sequence of domain 
losses Sm — jSq2~"^, and the error is small enough, the process can be iterated 
infinitely often and converges to a solution of the equation. The argument also 
shows that the hyperbolic splitting converges. 

Proposition 5.1. Choose an initial approximation Xq = 0, Kq, where Kq G ND{pq). 
Assume that Ko{Dpg), the range of Kq is at a distance r > from complement of 
the domain of definition of F . 

Assume {Xm-iT Km^i) G N D{pm-~i) is an approximate solution of equation ^ 
and that the following holds 

(84) ||i^,„_i-ifo||p,„_i <r/2, 

where r is chosen sufficiently small so that we can apply Lemma \4-9[ the constants 
in ND{pm~i) are chosen uniformly and that the range of Km-i is inside the do- 
main of definition of F. Assume furthermore that (|24p holds so that we can apply 
Lemma \4-^ 

Denote by C expressions that depend only on u, I, \F\c'i-{Br)j \\DKm-i\\p^_i, 
ll-'^if'^'" 1(9) llpm-i ' l(a'^s('3m-i))~"'"l C'''^d\(sMg{Am-i))~'^\ and, hence, can he chosen 
uniformly if Km-i is in a sufficiently small neighborhood of Kq as indicated in (j84p . 

Let Am-i, A,„__i be the corrections produced in Lemma \4-S\ 

IfEm—i is small enough such that 

(85) CKS-l\~'\\Em-i\\p„^_,<r/2 

then, the set (Km-i + ^m-l){Dp^_-^-s„^-l) ^s well inside the domain of definition 
of F and Em{d) = Tuj{\m,Km){0) satisfies (defining pm = Pm-i - "i^m-i) 

(86) \\Emh^<C^'5;^%\\Em-i\\l_,. 



Proof Wehave A„._i(0) = n^^_^(g)A,„_i(0)+n^^ _^(g)A,„_i((?), where n^^_^(g) 
is the projection on the hyperbolic subspace. Proposition l4.19l and Proposition l4.20l 
respectively, particularized to 5m-i give us that 

||A,„_i||p,„ <C«;2C-ill^^™-i||p„_, 
and using Cauchy inequalities. 

\\D/^m-i\\p^ < CK'd;;l\-'\\Em-i\\p„^_, 

Using ([55)) . and the previous estimates on Am-i, we see that the range of Km = 
Km~i + A„j_i is well inside the domain of definition of F so that we can define 
FoKm. 
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Define the remainder of the Taylor expansion 

n{X,X',K,K') ^T^{\K)-T^{\',K') 

-Dx,K:FU\K)iX-X',K-K'). 

Then we have 

Ernie) = +i?A,K^c.(A™-l,if„-l(0))(A,„_i, A„_i(0)) 

+ 7?.(Am-l, Am, Krn-1, Krn)iO). 

Using estimate (|75l) . for the error in solving the center equation and recalling that 
the equations on the hyperbolic subspace are exactly solved, we have 

Estimate (|86p then follows from Taylor's remainder bound 

|F o + A™„i)(0) - F o - DF o if„_i(0)A™_i(^)| 

<C\\D'F\U^,n^rm^. 

Note that, since ^„ go to zero, we can assume that the estimates from the Taylor 
remainder are larger than those from the error of the solution. □ 

5.2. Change of the hyperbolicity and the non-degeneracy conditions in 
the iterative step. The main goal of this section is to estimate the change of the 
non-degeneracy conditions in terms of the size of the error at the beginning of the 
iterative step. 

We begin by estimating the change in the invariant splitting. Later, we will 
estimate the change in the twist conditions. 

The first result Proposition 15 . 21 is a standard result in the theory of normally hy- 
perbolic sets that allows us to conclude that if we are given an approximately invari- 
ant splitting, which has some hyperbolicity, then there is a truly invariant splitting 
nearby. The proof is a reformulation in an a 'posteriori format of standard argu- 
ments on the stability of hyperbolic splittings [S574l IHP701 IFen72[ [PS991 [Pes04] . 
Since this will be part of an iterative procedure, we also need to obtain rather 
detailed estimates. 

As a corollary, we will obtain that, when we change the embeddings K in the 
iterative step, the change of the invariant subspaces will be controlled by the change 
in the embedding. Of course, since the twist conditions are just properties of the 
restriction of the derivative to an appropriate subspace, we will obtain that the size 
of the change in the twist conditions is controlled by the size of the change of the 
embedding. 

Notice also that Proposition 15 . 21 provides a way to verify the hyperbolicity out of 
a finite calculation and in particular, out of the results of a numerical calculation. 
We have also used Proposition [521 to identify the center space in Section 11.2.21 

Proposition 5.2. Assume that there is an analytic splitting 

(87) TK(e)M = © ® £^{9) 

which is approximately invariant under the co-cycle DF o K over T^^. That is, 
dist,pFoif(0)f^^," f^^-" )<5, 



34 



E. FONTICH, R. DE LA LEAVE, AND Y. SIRE 



where distp stands for the supremum of the distance when 9 belongs to Dp, the 
complex extension of the torus defined in We denote by H''''^'" the projections 
corresponding to the above splitting. 

Assume, moreover that, for some N ^ N, < fii, fl2 < 1, and some 1 < /ts, such 
that max(/ii,/i2) ■ jlz < I, we have 

\DFoK oT^~\e) X ••■ xDFoK{e)v\ <fli\v\ 
\DF-^ oKoT-^'^-^\9) X ••• X DF-^ o K{9)v\ <il^\v\ 



(89) 
and 



\DF oKo T^^^{e) X ■■■ X DFo K{9)v\ < fi^ \v\ 
(90) \DF-^ oKo rj(^"i'(6l) X • • • X DF-^ o K{9)v\ < fL^\v\ 

Assume that 6 < So, where So is an expression depending on N, \\DF o K\\p, 
\\DF-^oK\\p, ||n='^'«||p. 

Then, there is an analytic splitting 

TK(e)M = © © £^^(9) 

invariant under the co- cycle DF o K over T^^, which satisfies the characterization 
of hyperbolic splittings (jlip . p2p . (jl3p . 

The splitting above is unique among the splittings in a neighborhood of the orig- 
inal splitting of size Sq measured in distp. 

Furthermore, we have that 

j-gj^-j distp(5^"gj,f^"^P < CS, 

Imi,2,3 - Ml, 2, si < C'S, 
where C depends on the same quantities as Sq does. 

The previous result is applicable to all co-cycles over T^. It is important that 
the base is a rotation. As it is well known in the general theory of hyperbolic 
systems, if the base of the co-cycle had non-zero Lyapunov exponents, we expect 
that the invariant splittings are only finitely differentiable and not analytic even if 
the co-cycle and the base map are analytic. Some explicit examples are available 
in [dlLOlbj . 

In the statement of Proposition l5.2[ for typographical simplicity, we are assuming 
that the phase space is an Euclidean manifold so that we can compute the product 
DF o K{9 + uj)DF o K{9) and consider DF o K as a, co-cycle over T^. In case 
that the phase space is not an Euclidean manifold, the co-cycle is Sj^ 



where S is the connector introduced in Definition 13.91 This can be done provided 
that dist(F oi^(6'), K{9 + uj)) is small enough so that the connectors can be defined. 
The proof of Proposition 15.21 does not require any change beyond that to work in 
non-Euclidean manifolds. 
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Note that Proposition 15.21 implies immediately the persistence of invariant bun- 
dles under perturbations of the co-cycle. Given a co-cycle, its invariant bundles are 
approximately invariant under the perturbed co-cycle. The approximately invariant 
co-cycles can be obtained in many different ways, for example through numerical 
computations or through formal expansions. For the numerical applications we 
refer to HdlL06W . We also mention that [Mas05 ] computes Lindstedt series ex- 
pansions for quasi-periodic solutions in center manifolds for problems in celestial 
mechanics. These solutions are whiskered solutions in the full space and can be 
validated applying the results of this paper. 

Remark 5.3. Notice that the statements of the hyperbolicity conditions in Propo- 
sition 15.21 do not involve any constant Ct as in (fTTj) , (fT2l) , (fT3)) , but on the other 
hand, we include an N . From the point of view of mathematical theorems, both 
formulations are equivalent if we consider (jlip. (fT2| . (fT3| for fixed n = N . Note 
that \i jl > n and N are such that Ch{lJ-/p)^ < 1, the conditions in pill imply 
those in Proposition [521 The converse is trivial. 

We note that the constants Ch depend on the norm used in the space. Indeed, 
in theoretical applications, it is convenient to choose a norm such that (7^ = 1. 
Equivalently, one can choose a norm such that iV = 1. This indeed simplifies the 
notation. We have chosen not to take advantage of this simplification since the 
adapted norm is not commonly used in numerical applications. 

Remark 5.4. Another application of Proposition l5.2l that we will not develop here, 
is a bootstrap of regularity. If an invariant splitting is continuous, smoothing it, we 
obtain an approximately invariant analytic one and, applying Proposition 15.21 we 
obtain an analytic invariant splitting which has to coincide with the original one. 
See pohSO. ,HdlL06c] . 

Remark 5.5. With a view to the applications in [FdlLSOS] . we note that the 
arguments in the proof of Proposition 15.21 are rather soft (contraction mapping 
principle and such). Hence, they go through without changes when the bundles are 
Banach bundles. 

Remark 5.6. Notice that the proof of the existence of invariant subbundles given 
the approximately invariant ones is done one subbundle at a time. Hence, if we 
have two invariant subbundles (this is the situation considered in Proposition l4.16")) . 
the argument in the proof of Proposition 15.21 above leaves unchanged the invariant 
subspaces. Hence, the hyperbolicity constants /zi,/X2,/X3 and Ch in these spaces 
are unaltered. On the other hand, the projections on the invariant subspaces are 
altered because the projections depend on the splitting. The change of one of the 
subbundles changes all the projections. Of course, the change of the projections 
can be estimated by the change of the spaces, which is in turn estimated by the 
error in the invariance equation. 

The main application of Proposition 15.21 in this paper is the following result, 
Proposition 15.71 which estimates the change in the hyperbolicity hypotheses in an 
iterative step. 

Proposition 5.7. Assume that \\K — K\\p is small enough and hypotheses of Propo- 
sition [57l\ apply. Then there exists an analytic invariant splitting for DF o K. 
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Furthermore, there exists a constant C > such that we have the estimates 

(92) \\Ul'^^--U%^^^^\l<C\\K-K\\„ 

(93) \^l,~^i,\<C\\K-K\\p, z = l,2,3, 

(94) Ch = Ch. 

Proof of Proposition \5.SX The proof we present is very similar to the proof in |HdlL06b] . 
The ideas are very similar to the standard proof of the persistence of invariant 
splittings in [HPS77[ lPS99l IPes04| but we present them in an a-posteriori format, 
obtaining very quantitative estimates and we take advantage of the fact that the 
motion in the base is a rotation. This requires only some minor rearrangements of 
the argument in the above references. 
We will denote 

(95) ^Kio^-Jko) ^ 

^K{ff) — ^K(ff) ® ^K(0)- 

We clearly have 

(96) TK(e)M ®£l^,) 

and the splitting (j96p is almost invariant under DF o K. 

We consider the matrix of DF{K{9)) with respect to the splitting (|96|) : 

^nm) = (-;<^! -<^!^ 

The almost invariance of the splitting implies that ||ai2||p, ||a2i||p < Crj. 
We will construct the invariant subspaces corresponding to this splitting as 

graphs of linear functions u^{9) : Eg ^ £g and u'^{0) : Sg ^Re- 
computing the image of the point {x,u^{9)x), (resp. {u'^ {0)y , y)) and imposing 

that the images are in the graph of u^{9 + ui) (resp. u^{9 + uj)), we obtain that the 

graphs of m^, are invariant if and only if u^, v? satisfy 

(97) u^{9 + uj){a^x{9) + ai2(0)u'(e)) - a^^{9) + 022(^)^1(0), 

(98) an (0)^2(0) + a^^{9) ^ u\9 + cj) {021 {9)u^ {9) + 022(0)). 

As can be seen by elementary algebraic manipulations, equations ([97|) and ([98]) 
are equivalent to 

(99) u\9) ^ a^^\9){u\9 + uj){an{9) + ai2{9)u\9)) - a2i{9)), 

(100) u\9 + Lu) = {aii{9)u\9) + ai2(0))(a22(0) + a2i{9)u^i9))-\ 

We see that u^,u'^ are fixed points of the operators T^,T^ which are defined 
as the right-hand side of equation ((99|) and the rigth-hand side of equation (|100[) 
shifted by — w, respectively: 

T'[u']{9) =a22\0){u\9 + cu){an{9) + ai2{0y{9)) - a2i{9)), 

T^[u^]{9) ={aii{9 ~ Lu)u^{9 - to) + ai2{9 - Lu))x 

X (022(0 - w) + a2i(0 - cj)u^{0 - uj))-\ 

Now we concentrate on the operator T^. We introduce the space S — A{Dp, Ci) 
of analytic sections from Dp to the unit bundle of linear operators from f into 
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^K(e)^ i.e. the space of analytic maps u such that u{9) : £}(^g) — > ^^(g) is Hnear and 
< 1- Endowed with \\u\\s — sup^g^,^ ||u(6')||, 5 is a Banach space. Moreover 
S satisfies Banach algebra properties under the natural multiplications. 

We note that if ij is small enough and consequently ||ai2||, ||a2i|| are small, a 
reasonable linear approximation of is (obtained by eliminating all the terms 
that contain ai2, 0.21) 

[u^] {6) := 022' (9 + uj)aii{e) . 
An elementary computation gives 

= a^^iO) ■ ■ ■ 022^(6* +iN ~ l)uj)u\e + Nuj)aii{e + [N - l)uj) ■ ■■aii{e). 

Using the fact that is a quadratic polynomial operator, by performing alge- 
braic manipulations we obtain 

(101) max \\{TT[u'] iToT[u']\\p < Cry 
l|f llp<'? 

and 

(102) Ups^{iTY-(.%T)<Cv, 

where C depends on N. 

Note that, by assumptions ([55)) . ([90)) we have that if 77 is small enough, 

(T^)^ maps S into S. 

Then using (ISO)) . (I5D1) together with the previous estimates we have that 

(Tq)^ is a contraction from iS to S. 

This implies that also (T^)^ is a contraction. It is well known that then has 
a unique fixed point u in 5. 

Moreover the analyticity in G Dp is inherited by the fixed point of the contrac- 
tion. Hence u depends analytically in 9. 

Furthermore, we have the standard fixed point estimate 

(103) ||„||^<^((ri)^^(0)-0)<C77, 

1 — a 

where a — C"6 and the constant C depends on N,Ch, ^J■l, ^J■3■ This estimate gives 
that c?p(£^^^gp f^jgj) < C77. Then, since the spaces are Cjy-close we also 

have l/ii — /Lti| < C5. 

The proof so far, gives us the existence of invariant spaces as in ([95)1 . This clearly 
gives us the existence of the invariant bundle and the invariant bundle £'^(^gy 

We remark that exactly the same proof works if we take the splitting 

£K(e)^ £K(0)®£K(e) 
'-■KiyQ) — '-■K{ey 

Hence, we also obtain the existence of the bundles £^^1^^-^ and £\^Qy The invariant 
bundle f^^g^ is obtained as n £'^(Qy 



This concludes the proof of Proposition [ 

□ 

Proof of Proposition \ 5. 7[ We just observe that we can take the invariant splittings 
for DF o K a,s approximately invariant for DF o K. Using Cauchy estimates, we 
see that we can take 5 = C\\K — K\\p. Therefore, (|92| follows from estimating the 



38 



E. FONTICH, R. DE LA LEAVE, AND Y. SIRE 



change of the spaces. The conclusions (|93| . (|94l) follow from the observations in 
Remark 15.31 

□ 

The next Lemma |5 . 81 provides the perturbation for the remaining non-degenerate 
conditions. The idea is very simple. The twist condition is just the norm of a 
matrix obtained by restricting the derivative to the tangent and projecting it on 
the symplectic conjugate directions to the tangent. Cauchy estimates allows us 
to estimate easily the changes of these spaces. The estimate of the change of the 
derivative when we change the embedding is just the mean value theorem. 

Lemma 5.8. Assume that the hypotheses of Provosition [5A\ hold. // ||i?m-i ||p,„_i 
is small enough, then 

• If DK^-^_^DKm-i is invertible with inverse Nm-i then DK^DKm is in- 
vertible with inverse Nm and we have 

WNraWp^ < ||iV™_i||p,„_, +a„-l«'<5rn*_'r'^||i?™-l||p,„_,. 

• If a,vg{A,n-i) is non-singular then also &vg{A„i) is and we have the esti- 
mate 

|(avg(A„))-i| < |(avg(A„_i))-i| + C:„_i^25-^_\^+^)||i?,„_i||p„_,. 

• // avg ((5,„_i) is non-singular then also £tvg{Qm) is and we have the esti- 
mate 

|(avg(Q„0)-i| < |(avg(Q™_i))-i| 

Proof. For the first, we refer the reader to ( [dlLGJVOSj . Section 5) since the proof 
is identical. We turn to the second and third points. The estimates just come from 
writing Km ~ Km-i + using estimates (I25|) - (I26I) and neglecting quadratic 

error terms at the price of changing the constants. □ 

5.3. Convergence of the scheme. It is by now classical that, under sufficiently 
strong smallness assumptions, the iterative scheme can be iterated indefinitely and 
that it converges. Similar arguments can be found in almost any paper in KAM 
theory, in particular |Zeh75j . |Zeh76| . [Mos66bj . |Mos66a| . [dlLOla] . The notation 
in this paper matches closely that in [dlLGJVOS] so that the modifications, at this 
stage are rather minimal. 

Recall that we have identified a set of embeddings in which we can obtain uniform 
constants in the Newton step, see Proposition 15. II 

In the following Lemma 15.91 we show that, with the choice of domain losses 
given in (|105p if the initial error is small enough, the iterations do not leave the 
neighborhood where we have uniform estimates and converge to a solution of the 
problem, which also has hyperbolic splittings. 

Lemma 5.9. Using the previous notations, let Cm be the sequence of positive num- 
bers defined above. For a fixed < (5o < min(l,po/12) define for m > 0, 

(105) 5m = 5o2-", 

Denote pm = Pm-i - G5m-i and em = ||-Em||p,„. 

There exists a constant C depending on I, v, \F\c2(^g^-^, \J\c'^(B^)j ll-C^ollpo; 
ll^ollpo, |(avg(Qo))"^|, |(avg(ylo))"^|, l|n^n'fe)llpo; ll^llpo such that if the error 
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satisfies the following inequalities 

C2^''K'^So'^''eo < 1/2 

and 

then the modified Newton step can be iterated indefinitely and we obtain that Km 
converges to a map K^c G Ap^-QS^ which satisfies the non- degeneracy conditions, 
in particular, it is hyperbolic, and 

Moreover, there exists a constant D > depending on I, v, \F\ij2.(^Br)j \J\c^(Bt.)} 
\\DK^\\p„ IITVoIIpo, |(avg(Qo))-^|, |(avg (Ao))-i|, lln^^;^",) and ||G||p„ such that 

\\Koo - -?fo||po-6<5o < Dk^5q'^''\\Eo\\po- 

Remark 5.10. We note again that by the vanishing lemma 14.91 the sequence 
{'^n}n>o converges to as n goes to +oo. 

Proof. As mentioned in the introduction of the section, the argument is quite stan- 
dard. 

To ensure that we can perform steps with the estimates in Proposition 15.11 we 
just need to verify that we do not leave the neighborhood of A'o given by ([84]) and 
that we satisfy the bounds ((55)) . 

We note that, in a concise notation. Proposition 15. II leads to the bounds 

With the choice 6m = Sa2^"^, we see that if we can perform m steps, we have: 

< ((7K4^-4iy^l+2+---+2"-i24'^[(™-l)+2("i-2)H h2'""^]g2'" 

^ ('/^o4i' „4 r-4iy, \2™ 

< [C2 K dg eo) , 
for m > 1, where we have used that 

(m-l)+2(m-2)+- • •+2"-2 = 2'"-^[{m^l)2-^"'-^^ + {m~2)2-^"'-^^+- ■ ■+!] < 2"^. 

We see that if eq is small enough, then, emi^"^" is so small than the conditions 
([85]) are true for the next step. Indeed, we note that the smallness conditions that 
we need to impose in eo are independent of m. 

Furthermore, we also observe that we also have Km — Kq ~ X]"io^ ^i- Hence 

m — 1 m — 1 

\\Km~Ko\\p,^^ < ||A,||,. < CK'{CKHo"'^ofSo'-2^'-'. 

1 = i=0 

We note that, by taking eg small enough we can make the right-hand side of the 
last formula as small as desired uniformly in m. In particular, by taking eo small 
enough we can ensure the assumption ([84)1 for all m. 

Therefore, if we assume that eo small enough, we can ensure that we can repeat 
the iterative step infinitely often and that the iteration never leaves the neighbor- 
hood identified in (l84l). 



40 



E. FONTICH, R. DE LA LEAVE, AND Y. SIRE 



We also note that we have 



oo oo oo 




oo 



< Dk- 







The absolute convergence of the above series shows that Km converge to a limit 
and the last bound establishes the conclusion (|19p. We note that since we had 
assumed (j84p , we have that Koo admits a hyperbolic splitting. Since the change in 
the hyperbolic splittings is bounded by the change of the embedding (see Proposi- 
tion l5.7p . we see that the hyperbolic splittings also converge to the limiting one. □ 



In this section, we prove Theorem l3.14l We closely follow the proof in [dlLG JV05] . 
Similar results are more or less implicit in the treatment of whiskered tori in 'Zeh76] . 
For fully dimensional tori local uniqueness results appear in [Mos66b, SZ89, Sal04j. 
As we have argued before, local uniqueness results allow us to deduce results for 
flows from results for maps. 

The proof of Theorem 13.141 is based on showing that the operator DT^ (K) has 
an approximate left inverse (as in |Zeh75] ). Notice first that the composition on 
the right by every translation of a solution of ^ is also a solution. Therefore, 
one cannot expect a general uniqueness result. Moreover, the second statement 
in Lemma [4.21 and the calculation on the hyperbolic directions show that, roughly 
speaking, two solutions of the linearized equation differ by their average. Moreover 
this difference is in the direction of the tangent space of the torus. 

The idea behind the local uniqueness result is to prove that one can transfer the 
difference of the averages of two solutions to a difference of phase between the two 
solutions. 

Now we assume that the embeddings Ki and K2 satisfy the hypotheses in Theo- 
rem l3.14l in particular Ki and K2 are solutions of ([6]), or ([7]) with A = 0. If r 7^ we 
write Ki for KioTt which is also a solution. Therefore !Fui{Q, Ki) = .7-^(0, K2) = 0. 
By Taylor's theorem we can write 



6. Proof of the local uniqueness theorem 



(106) 



= K,) - K2) =Dx,k:FUO, K2){0, - K2) 

+ n{0,0,Ki,K2), 



where 



TZ{0, 0, Ki,K2)^^ [ D^F{K2 + t{K^ - K2)){K, - K2f dt 



Then, there exists C > such that 



\\n{0,0,K,,K2)\\p<C\\K,-K2\\l. 



Hence we end up with the following linearized equation 

(107) Dx,kTUO, K2m Ki - K2) = -7^(0, 0, Ki,K2). 



We denote A = i^i - 
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Projecting ()107|) on the center subspace with ^'k2{8+uj)^ writing A'^{6) — U'j^^^g~^A{9) 
and making the change of function A'^(6') = M{d)W{9), where M is defined in 
with K = K2, we obtain 

DF{K2{0))M{9)W{e) - M{9 + uj)W{9 + uj) 

(108) =-n^^(«+^)7e(o,o,Xi,X2)(e). 

We note that since K2 is an exact solution, range M (9) coincides with S'^^^gy See 
Subsection 14.1.31 

Applying the property DF{K2{9))M{9) = M{9 + uj)S{9) for solutions of jS]), 
multiplying both sides by [M^ J{K2)\{9 + uo) and using that J{K2)M is invert- 
ible we get 

S{9)W{9)~W{9 + u) 

= -[{M^ J{K2)Mr^M^ J{K2)]{9 + c.)n^^(,+^)7e(0, 0, K^, K2){9). 

Since W solves the previous equation, we get bounds for it using the methods in 
Section r4. 2. 41 We write W = (Wi, W2). Since S is triangular we begin by looking 
for W2- We search it in the form W2 — + ayg{W2)- We have ||W2-^||p_5 < 
CkS-''\\Ki - Xallp- For Wi we have 

W,{9)-Wi{9 + uj) -T2(^)(^^^(,+^)7^(0,0,i^l,i^2))l(^) 

(109) - Ai9)W^{9) - Ai9)aYg{W2), 

where T2 = Nj DKj J{K2y^[DK2N2DKj - Id] J(/-i:2) and A^2 = DKJDK2. 

The condition the right-hand side of (|109p to have zero average gives |avg (1^2) | < 
CkS-^WKi- K2\\l. Then 

\\Wi - avg(l¥i)||p_25 < Ck'^S-^'^WKi -K2\\l 
but avg{Wi) is free. Then 

||A^ - (avg(A^)i,0)T||p_25 < Ck'S-'-'WK, - K2\\l 

The next step is done in the same way as in }dlLG JV05] . We quote Lemma 14 
of that reference using our notation. It is basically an application of the standard 
implicit function theorem. 

Lemma 6.1. There exists a constant C such that if C\\Ki — K2\\p < 1 then there 
exists an initial phase ti E {t eMJ \ \t\ < \\Ki — K2\\p} such that 

avg(r2(0)n^,(e)(i^i °T., -i^2)(e)) = 0. 

The proof is based on an application of implicit function theorem in R'. 
As a consequence of Lemma |6.1[ if ti is as in the statement, then K o T^-^ is a 
solution of ([6]) such that if 

I¥ = [Af^ J(if2)A7](^? + cc;)-i[M^^(^2)](0 + c^)n^^(e)(i^i o T,, - ifa), 
for all 5 E (0,p/2) and we have the estimate 

\\w\\p-25 < CK's~^''\\n\\l < CnH-^''\\Ki - K2\\l. 

This leads to on the center subspace 

mhie)iKioTr, - K2)\\p^2S < CnH-^^\\K, - i^2ll^. 
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Furthermore, as in Section 14.31 taking projections on the hyperbolic subspace, we 
have that A'' — U'^^i^g-^{Ki — K2) satisfies the estimate 

\\/:^%-25<c\\np- 

All in all, we have proven the estimate for Ki o T^-^ — K2 (up to a change in the 
original constants) 

lli^l O Tr, - K2\\p-2S < Ck^S'^^WKi ~ K2\\l. 

We are now in position to carry out an argument very similar to the one used in 
Section [Ol We can take a sequence {Tm},„>i such that |ri| < Iji^i — K2\\p and 

km - Tm-l| < 11^1 O 7V„_i - K2\\p^_i, m > 2, 

and 

lli^i o Tr„^ - K2\\p„^ < CnH;^^^\\K^ O Tr^_^ - K2\\l_^, 

where i5i = p/4, 6m+i = ^m/2 for m > 1 and po = p, Pm = Po- ^2^=1 for m > 1. 
By an induction argument we end up with 

\\Ki o - K2\\p^ < {CkH^^-2^''\\K, - XsllpJ^-a-^-^™. 

Therefore, under the smallness assumptions on — ii'2||po > the sequence {Tm\rn>i 
converges and one gets 

\\K,oTr^-K2\\p/2=0. 

Since both Ki o T^^ and K2 are analytic in Dp and coincide in Dp/2 we obtain the 
result. 

7. Applications 

In this section, we collect several consequences of our main theorem. We note 
that these consequences follow mainly from the fact that we have formulated the 
theorem in a posteriori style without reference to an integrable system. 

7.1. Lipschitz dependence vifith respect to the frequency. Estimates of the 
measure occupied by the tori. The basic idea is that if we have an embedding 
K that solves the equation for one frequency, then it solves approximately the 
equation for a nearby frequency. Then, applying Theorem 13.111 there should be a 
solution for a new frequency which is close to the the original one. Performing the 
argument with care, we see that this implies Lipschitz dependence of the solution 
on the frequency. Similar ideas were indicated in [Zeh75j . We remark that this 
Lipschitz dependence leads to estimates on the measure occupied by the tori in the 
perturbative case. We concentrate on the case of maps since, as we have shown, it 
implies the corresponding result for flows. 

We assume that F is defined and analytic in a sufhciently big complex domain of 
an Euclidean manifold Ai. We consider u) € D{k,v) with k and v fixed and we sup- 
pose that G Ap satisfies ^ and is non-degenerate in the sense of Definition [3111 
For all ijj' G D{k, v) we have 

and therefore, applying the mean value theorem and Cauchy estimates, we have 

(110) \\FoK^-K^oT^,\\p_s < CS-^\\K^\\p\lu-lu'\ 

for all 6 £ (0,p). If jw — a;'| is small enough, namely Ck^S^^''~^\\Kuj\\p\uj — uj'\ < 1, 
applying Theorem l3.11l we obtain that there is an embedding Koj' satisfying ([6]) with 
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the frequency w'. Furthermore, taking the value of 5 in Theorem lS. in appropriately, 
one gets 

\\K^ - K^'\\(p-5)/2 < Ck\p ~ Sy^^S-'WKJplu; - u;'\. 

We note that, in the domain of applicability of the previous argument, the Lip- 
schitz constant is uniform since we are assuming that k and v are fixed. Since 
the set of uniformly Diophantine vectors D(k, v) is locally compact, we can cover 
a bounded subset of D{k,, v) with a finite number of balls in which the previous 
argument applies and therefore in this set we get Lipschitz dependence on lo. 

Moreover, the frequency of is given by the formula (up to some multiplicative 
constant) 

(ill) uj = Tl^ j {Fok^{e)-k^{6))d9, 

where F and K are the lifts of F and K to the universal cover of M and 11;^ is the 
projection over the lift of the angle variables. 

Thanks to formula (|llip . it is straightforward to see that the map i-^ uj 
is Lipschitz. Hence, we conclude that the map w i— > is bi-Lipschitz from the 
set of Diophantine vectors with fixed Diophantine constants. Since the set D[k, v) 
has positive Z-dimensional measure, we conclude that the set of tori also has 21- 
dimensional measure, i.e. 

n'-'i y if^(T'))>o, 

where 7i^' stands for the Hausdorff measure. 

7.2. Analyticity with respect to parameters. The proof of the existence of tori 
associated to a fixed frequency lo presented here leads to analyticity in the depen- 
dence with respect to parameters. Later, we will see that this leads to analyticity 
properties of some series expansions, such as Lindstedt series. The argument is 
already contained in |Mos67j . However the argument presented here is somewhat 
simpler than the one presented in that reference. 

Given a family of functions F^ and a family of approximate solutions both 
depending analytically on parameters G C/ C (p > 1) and continuous in [/, 
we see that the assumptions of Theorem 13.111 are satisfied uniformly in G C/. 
Consequently, there is a true solution nearby which also depends analytically on 
the parameters 77. 

The proof is very simple; we just observe that the iterative step is analytic for 
77 6 J7 (resp. continuous for 77 S C/) if the family and the error are. 

The procedure and the estimates for one step of the iterative procedure are 
stated in Lemma 14.21 and in a more detailed way in the statements and proofs of 
Proposition 14. 191 and Proposition 14.201 We just call attention to the fact that the 
correction applied at each step of the application of Proposition l4. 191 relies on some 
explicit algebraic formulas — involving derivatives with respect to 9 — and to use 
the solution of some small divisor equations. Note that the solution of the small 
divisor equations is obtained applying a linear operator which is independent of 
77. Also the method of obtaining the projection on the hyperbolic directions done 
in Section 14.31 and summarized in Proposition 14.201 preserves the analyticity on 
parameters since A*'" are obtained as sums of uniformly convergent series. 
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Clearly, the analyticity properties with respect to rj are preserved by all these 
steps. Hence, the corrections applied in one step of the iterative Lemma l4 . 2 1 depend 
analytically on parameters when the error does. Of course, the error depends 
analytically on 77 e [/ (resp. continuously on 77 e [/) if the approximate solution at 
the start of Lemma |4. 1 71 depends analytically on the parameters, since to compute 
the error from the approximate solution, we just have to compose with the function 
Fri and translate. Therefore, we conclude that the application of Proposition [47191 
preserves the analyticity properties with respect to parameters. 

Hence, in all steps of the iterative process used in the proof of Theorem l3.111 the 
functions {^m}„igN depend analytically on parameters ranging on the open set U 
(and continuously on the boundary). Of course, in the iterative step, we decrease 
the analyticity domain in the variables 9, but not the analyticity domain in rj. 

We also observe that in the proof of Proposition I4.19[ the estimates on the 
correction applied at each step depend only on the sizes of the error and the non- 
degeneracy conditions. Also, we observe that the estimates on the change of the 
non-degeneracy conditions are uniform on the size of the corrections. In particular, 
if we assume that the smallness and non-degeneracy conditions hold uniformly for 
?7 € ?7, we can apply Lemma 15.91 to obtain that there is a sequence of analytic 
functions in 9, rj converging uniformly for 9 G Dp^ and r] £ U. 

In the paper [dlLOOOj there is an alternative point of view for results with pa- 
rameters. One can apply an abstract KAM implicit function theorem as in j Zeh75] 
to spaces of analytic functions in other Banach spaces. These kind of arguments 
were used to deal with rather degenerate problems. A more detailed study of KAM 
theorems with parameters appears in |Van02j . 

7.3. Small twist theorems and small hyperbolicity theorems. Small twist 
theorems have been introduced in |Mos62[[Kyn68| to deal with problems in celestial 
mechanics. The idea of small twist (and small hyperbolicity) theorems is to give 
conditions that ensure the convergence of the Newton-type method even if the twist 
is close to be degenerate. It goes through a more precise analysis of the constants 
involved in the Newton scheme. 

We refer to [Mos62[ [KynSSI IOrt97[ [Ort99| for applications of small twist theorems 
to celestial mechanics and to the stability of oscillators. 

The goal of this section is to provide, as a corollary of the proof of our main 
Theorem 13. Ill a small twist and small hyperbolicity result. 

By examining carefully the proof involved in the iterative step in the KAM 
method (see Section |4. 2. 4|) . we get the following proposition. 

Proposition 7.1. • There exist two positive numbers a, (3 such that the con- 

stant C in equation (j75p (which here will be denoted C) depending on I, 
K, V, |F|c2(B.), \\DK\\p, \\N\\p, |(avg(A))-i|, |(avg (Q))-i|, HH^^.^t and 
||G||p can be estimated by 

C^<||H^^(,)||pmax(l,pX||p)"max(l,l|A^l|,)^(l(avg(A))-i| 

+ |(avg(Q))-i|), 

where A, Q are defined in (I16|) . psp respectively. 
• The constant C in equation (|8ip (which here will be denoted C^) depending 
on the hyperbolicity constant /ii (resp. ^2), the norm of the projection 
||n^,gJ|p (resp. ||n^,gJ|pJ and \\G\\p and the constant Ch involved in J77]) 
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and can he estimated by 

(113) < a(l + C=)max(||n^(,)||p-^, ||n^(,) ||,-^). 

• As a consequence of the two above items, the constant C appearing in The- 
orem \S.11\ ( the one which enters in equation ) can be bounded by 

(114) 

c>ax(||n^(,)||p-^,||n^(,)||p-^)2 

+ ||nj,^(,)||>ax(l, \\DK\\,r max(l, \\N\\,f {\{^yg{A))-^\ + |(avg (Q))-i|)^ 

The argument presented in this paper gives a = 4,/3 = 2, but there are other 
variants of the argument which give better values. We have not optimized the 
bounds. 

To prove Proposition !?. II we note that to find A*'" we just use formula ([55]) from 
which the claim follows by estimating the sum using the triangle inequality and 
using the sum of the geometric series. 

The estimates claimed for the constants related to follow by observing that 
the solution is obtained by applying the following operations: multiplying by the 
matrices M, M^, multiplying by the matrices N , modifying the constants A and 
choosing the average of W2- The latter steps are estimated by multiplying by 
|(avg(Q))-i|andl(avg(A))-i|. 

We also recall that the remainder of the Newton method is estimated by the 
remainder of the Taylor expansion. Hence, it is estimated by the square of ||A||p_5. 

Let m > be an index for the Newton step and denote C the constant involved 
in the third item of Proposition 17.11 We have for some u > 

where p„ = /7„,_i - (5„i_i. 

It is standard in KAM theory (see Theorem 13. lip that if 

then the Newton method with 5m = 2~™5o converges to a solution. 

Therefore, even if the twist and the hyperbolicity are close to degenerate so that 
C is large, if the initial error is small enough, one gets convergence of the scheme. 

Small hyperbolicity arises naturally in perturbations of integrable systems. The 
integrable system, of course, has no hyperbolic behaviour, but an averaged sys- 
tem has some small hyperbolicity. Indeed, similar considerations for periodic orbit 
happen already in [Poi991 Ch. 74, 79]. 

The papers |JdlLZ99l ldlLW04[ IC W991 IChe99[ IeIiMI ITre94| consider perturba- 
tions of integrable systems at resonances, where the hyperbolicity is small and the 
present result can be applied. These papers differ in several important aspects such 
as the dimension and the topology of the tori considered. The methods are also 
different. 

All of the above papers consider perturbations of quasi-integrable systems i?o -I- 
sHi. 

The paper [JdlLZ99j shows that, given some appropriate non-degeneracy con- 
ditions on the perturbation, it is possible to construct formal series of approxi- 
mate solutions in powers of e. Truncating the series up to order N , it is shown 
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that the error of the power series can be bounded by CN°- e . Similarly, for 
some of the solutions, H^^g*^ are of order (Ree)~^/^, the hyperbolicity constants 

(1 — Aii)^^, (1 - A*2)^^ are of order (Ree)"^/^ but (avg {A))^^ and (avg (Q))^^ are 
of order 1. The Diophantine constants can be assumed to be fixed. 

If we fix r > sufficiently small and consider the set r < jej < 2r, we can choose 
the order of truncation so that the error is less that Cexp(— be""^). Then, the small 
hyperbolicity result applies to show that there are invariant tori, which depend 
analytically in e for e such that r < |e| < 2r and Ime > Cexp(— 6(Ree)'^). Since r 
is arbitrary, we obtain that there are hyperbolic invariant tori in a ball except, at 
most in a wedge around the positive real axis which is exponentially thin. 

We refer to |JdlLZ99) for precise conditions on the series so that we can get 
the perturbation series as above. For the case of two degrees of freedom, the 
paper [JdlLZ99] considers the existence of elliptic tori. We note that the method 
of |JdlLZ99j is based on reducibility. This paper shows that we do not need to 
use reducibility for the hyperbolic directions. The study of elliptic directions has 
experienced very significant progress in the last years, but we will not mention it 
here. 

The paper [dlLW04] considers also weakly hyperbolic tori around periodic orbits 
generated by resonances. Note that the tori considered in |dlLW04j are secondary 
tori, that is tori that cannot be deformed to tori with the same frequency in the 
unperturbed system. Indeed, the tori can be deformed into tori with less angles. 
Since [dlLW04| involves reduction to a center manifold it only concludes that the 
tori are finite differentiable even if the system is analytic. A result which is improved 
by using the results provided in this paper in Section [7.51 

7.4. Secondary tori and whiskered tori close to rank-1 resonances. The 

method described in this paper can accommodate to study secondary tori. Indeed, 
the development of algorithms which could deal with secondary tori was an im- 
portant motivation to modify the invariance equation by adding a term containing 
A. 

We recall that secondary KAM tori are invariant tori, such that the motion on 
them are conjugate to an irrational rotation but in contrast to the usual KAM tori 
which are homotopic to T' x {0 G R^''"'}, the secondary tori are homotopic to 
T'-*^ X {0 e R2d-;+fc|_ 

The existence of secondary KAM tori is very apparent in numerical explo- 
rations. For example, in two-dimensional maps, they are known as islands. In 
two-dimensional maps, islands are quite visible and they may occupy a large mea- 
sure of the phase space. 

Note that secondary tori are not present in the integrable system and their 
existence is not guaranteed by the standard pertubative KAM theory, which is 
concerned with the persistence of the invariant tori already present in the integrable 
system. In contrast, they are generated by the perturbation. The perturbation 
theory is somewhat unconventional since the unperturbed system does not present 
the phenomenon. Perturbative proofs of existence of secondary tori are done in 
[dlLW04j and in |DdlLS06|. 

In the recent papers [DdlLSOSl IDdlLSOH] it is shown that these secondary tori can 
be used as effective tools to generate diffusion and, in particular, to overcome the 
large gap problem in the study of diffusion. The paper [HdlLOO] argues heuristically 
and verifies numerically that, in multiparticle systems that will be considered in a 
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follow-up of this paper, in particular in the celebrated Fermi-Pasta-Ulam |FPU55j . 
the secondary tori occupy a much larger volume of phase space than the primary 
tori. 

The method to construct whiskered tori in [dlLW04j was to show that, under 
explicit conditions on the perturbation, the rational frequencies give rise to periodic 
orbits, some of which admit center manifolds. Under appropriate non-degeneracy 
conditions, these center manifolds contain tori which are invariant under the re- 
striction. These invariant tori, are whiskered tori for the full system. They are 
secondary since the directions corresponding to the center directions can be con- 
tracted to a point. 

The paper |DdlLS06] shows that secondary tori are generated by resonances in 
systems such that the unperturbed system has a two-dimensional normally hyper- 
bolic manifold. The method of proof is to show that, near the resonances, one 
can approximate the system by a system which is pendulum like. This pendulum 
has orbits that are rotating. In [DdlLS06] . it is shown that one can consider the 
real system as a perturbation of the pendulum and, therefore that some of the tori 
present in the pendulum are also present in the real system. See also |DH06| . 

One of the difficulties of the method in [ DdlLS06] is that the action variables 
near the separatrix are singular. This difficulty is, of course, not a problem for the 
method developed in the present paper. The method used in [DdlLS06] to overcome 
the singularity of the action variables was to perform more averaging steps, which 
required assuming more regularity of the perturbation. Applying the methods of 
this paper allows us to reduce the number of derivatives assumed in [DdlLS06] . 

As we will discuss in more detail in Section [7.51 by using reduction to center or 
normally hyperbolic invariant manifolds, one can only prove that the obtained tori 
are finitely differentiable even if the mapping is analytic. Using the results of this 
paper, we will show that these tori are actually analytic if the map is. 

7.5. Bootstrap of regularity of invariant tori. In this section, we show that if 
an analytic exact symplectic map F admits an invariant torus, with the maximal 
number of hyperbolic directions permitted by the symplectic structure, of class 
with r large enough, then the torus is actually analytic. Similar results for 
Lagrangian tori have been proved in [SZ89| . 

Proposition 7.2. Let F : A4 A4 be an analytic exact symplectic map. Let 
LO G -D(k, v) for some k > and v >l, and K ^ M satisfy 

(1) K is a solution of the equation F o K — K o = 0. 

(2) K is non- degenerate in the sense of Definition \3.4\ 




Then K is analytic. 

Remark 7.3. One case when Proposition l7.2l is useful is when the tori are produced 
by a reduction to a center manifold or a normally hyperbolic manifold. These in- 
variant manifolds are only finitely differentiable. Applying the above result. Propo- 
sition [7?2l we can conclude that the tori are analytic. 

The paper |dlLW04j constructs whiskered tori by applying the KAM theorem for 
Lagrangian tori to the restriction of the system to a center manifold. The papers 
[DdlLSOOt IDdlLSOSt IDdlLS06] consider tori in a normally hyperbolic manifold. In 




(3) K IS C with 



(115) 



r > Aiy. 
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particular, |DdlLS06] constructs secondary tori. We conclude that, in the case that 
the considered system is analytic. Proposition 17. 21 shows that the tori are analytic. 

The idea of the proof is very simple. We approximate the function K by an 
analytic one which will be an approximate solution of equation ([6]). Applying our 
main Theorem 13 . 1 II we will obtain that there is an analytic invariant torus nearby. 
The smoothness in the assumption enters because the hypotheses of Theorem l3.11l 
involve the size of the error in a complex strip. The uniqueness result Theorem l3.14l 
will give that the analytic torus obtained this way coincides with the original one 
up to a translation in the "angles" . 

The construction of the analytic approximations could be done in many different 
ways. For example, truncating the Fourier series of K would do, if one assumes a 
condition stronger than (|115p . 

As it is well-known since [Mos66bj , a very efficient way of approximating smooth 
functions by analytic ones is performing a convolution with a suitable kernel. 

Following [Zeh75|, IMos66b] , we introduce smoothing operators that provide nat- 
ural ways of approximating smooth functions by analytic ones. 

Definition 7.4. Let u : M} ^ R he a C°° even function identically 1 in a neigh- 
borhood of the origin and with support contained in the unit hall. Let it : K' — > C 
he the Fourier transform of u and denote by v the holomorphic continuation of u. 
For f e C"(R') andt>0 we define 

(116) St[f]iz):=t' [ v{t{y-z))f{y)dy. 

Jm' 

The map St defines a linear operator from C°(T') to the space of analytic maps 
from Dp to C, p > 0. Moreover, St is an analytic smoothing operator in the sense 
of [Mos66b| IZeh75| since it satisfies the following proposition (see |Zeh75| Lemma 
2.1] for a proof). We recall that if g G Ap, \\g\\p = sup^g^i^ |g(-z)|- 

Proposition 7.5. Let r G M+\N. There exists a constant ni — Ki{l,r) such that 
for all t > 1 and all f G C"^(T') we have 

(1) |(5*-Id)[/]|co </ti|/|c'- t-\ 

(2) \\St[f]\\t-i < Kilflco, 

(3) \\{Sr-Stm\\r-^ <ti,\f\c^t-^, forallT>t. 

We note that, since the smoothing operator commutes with derivatives, we also 
have the following extensions of H]) and ([2]) for s < r 

(117) |(^t-Id)[/]|c'- <Ki|/|c. ^-'■+^ 

(118) \\D^St[f]\\t-i < ^i\f\c^. 

Proof of Proposition \ 7. S\ We consider S** [i^T], the smoothed version of K with t>l. 
Our first goal is to estimate the error in a domain of size t~^£^ with ^ G (0, 1). 

We note that, by © in Proposition [731 and pT5)l . H^f < ki\K\co and 
ll-DS'f [XJIIt-i < Ki|_ftr|(^i remain bounded uniformly in t. 

Lemma 7.6. For t>\ and f G C"'(T') we have 

I \\D'St[f]\\t-i - \D'f\co I < 2n^\f\c^t-\ < s < r - 1. 

Proof. Since / G C''(T'), St[f] is analytic and T' and Z^j-i are compact, there 
exists xo G T' and zq G A-i such that IDVIc" = \D''f{xo)\ and \\D'' St[f]\\t~^ = 
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\D''St[f]izo)\. Assume that ||i:i'^S't[/]||t-i > \D'f\co. Then applying ^TJii and 
(fTTS]) we have 

0<\D^St[f]izo)\^\D^f{x,)\ 

< \D'St[f](zo) - D'St[f]{Rezo)\ + \D' St[f]{Re zo) - f{Rc zo)\ 

+ \D'f{Rczo)\-\D'f{xo)\ 

<^tl|/|c.+l^-l+Kl|/|c^<-''+^ 

If ||_D'*S't[/] llf-i < l-D^/Ico we argue in a symmetric way and we obtain the resuh. 

□ 

bmce X(T') is real we have that for t large enough, St[K]{Dt-i) is contained in 
the domain of F. Therefore, if t is large enough, we have that ||F o 5't[ii']||t-i and 
ll-F o S't[i4r] — S'([-ftr] o Ttj||(-i remain uniformly bounded. 

On the other hand, by ([T]) in Proposition 17.51 and the fact that K satisfies the 
functional equation we have that 

|F o St[K] - St[K] o T^lco <\FoSt[K]-Fo K\co + \St[K] oT^-Ko T^\co 

< Ki\K\c-{\F\ci+l)t-''. 

Therefore, using the interpolation inequality ([4]) in Proposition l2.5l with pi — t^^ 
and p2 = 0, we obtain that 

\\FoSt[K]~St[K]oTjt-,^ 

(119) < \FoSt[K]-St[K]°T^\c'o^\\FoSt[K]~St[K]oTjj_, 

Since all the non-degeneracy constants involve the first derivatives, by Lemma [7.6l 
we can perform the perturbative arguments in Section O 

The constants in the non-degeneracy assumptions remain uniformly bounded for 
St[K] in a neighborhood of size and, a fortiori, in a neighborhood of size t~^£^. 

Therefore, we can apply Theorem l3. Ill with po — t^^£, and 6 = t^'^£,/l2 provided 
that we can find t > 1 such that 

^(^-l^)-4^^-r-(l-C) ^ ^ 

for some constant C > 0, which depends on I, \\DSt[K]\\t^i^, 
(avg (yl))^^ I, |(avg (Q))^^|. By Lemma [7.61 if ^ is big enough, the constant C can 
be chosen independently on t. 

The condition r > Av implies that there exists ^ close to and t sufficiently 
large such that the previous inequality holds. Applying Theorem 13 . 1 1 1 with initial 
approximation Kq = St [K] we conclude that there exists an analytic solution 
of equation ([6]) defined on Di-i^ji which satisfies 

where Ci depends onl,v, \\DSt[K]\\t-i^, \\N\\t-^^, \\A\\t-i^, |(avg |(avg (Q))-i 

As before Ci can be taken independent on t. From (3) in Proposition [73] we have 
{r>t) 

\\Sr[K] - St[K]\\r-i^/2 < \\Sr[K] - St[K]\l-l < C^r'^' 

with C2 independent on t. 
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We will apply Theorem 13.141 with Ki and K2 being and respectively, 
with t, r > 1 . The application of this result requires Condition ([20]) which in our 
case reads 

(120) C3«:^(^^)-^n|i^r - K^r-^m < 1. 

The constant C3 depends on/, I/, \\K^\\,-i^/2 < ||i^rllt-i?/2, ll^t||t-i? Pt||t-i«, 
(avg(At))~^|, \ {avg{Qt))~^\, where Nt, At and Qt are the expressions introduced 
in Definition 13.41 corresponding to K"^ . As before C3 can be chosen independently 
on t,T E [1, 00), if t is big enough. We write C3 — S^'^k^C^. 

Lemma 7.7. There exists t > 1 such that if t > t there exists (ft.r G T' such that 

(121) K^oT^^^^K^. 

Proof. Using the previous notation we take t big enough such that the constants 
Ci, C2 and C3 are independent on t and such that 

(122) C32^T^''{2Cit^''-''^^-^'> + C2t^''-'') < 1. 

We define <,„ = 2"H, m > 0, and we claim that for tm < t < 2tm = t„i^i there 
exists (pt^.T <= T' such that 

Indeed, we apply Theorem 13 . 141 with Ki — and K2 = . We have 

\\KZ - K^\\r-H/2 <\\KZ - ^t„milt-e/2 + WStAK] - Sr[K]\\,-.^,2 

<C^{t-:0-V-^^ + C2t^: + Ci(r"iO-"^r-(i-«). 
Using that r < 2i„i, Condition (|120p is implied by 

which holds true by (|122p since i,„ > t. 

If r > t there exists fc > such that tk < t < tk+i- From the claim we can 
define (pt,T = J2,n2o ^t^,t.^+t + 'Ptfc.r- Clearly ipt^r satisfies (|121l) . □ 

Now consider Tj > < going to cxd. Since ft.Tj G T' there exists a convergent 
subsequence, which we denote again ft.r-: with limit (^00 G T'. 
Then 

- Xbo < - 5., [A1 bo + \Sr, [K] ~K\co 

Also, using that — o T^^^^. we get 

\Kr o T^^ - if Ico < |Xr ° T^oo ~ ° I^^*,., Ico + - K\co. 

Finally, taking limit as j goes to 00 we get K — oT^^ and hence K is analytic. 

□ 



7.6. Nontrivial stable and unstable bundles. 
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7.6.1. General comments and classification of bundles. In this section we describe 
some examples of whiskered invariant tori with non-trivial stable/unstable bundles. 
Theorem 13.111 applies to these tori while other methods in the literature do not 
seem to apply. We note that, for some systems (see |HdlL07| ). non-trivial bundles 
appear naturally when the systems experience resonances. We think that the study 
of bifurcations of the bundles of invariant tori deserves further exploration. 

We are very grateful to Prof. R. Gompf for very enlightening discussions and, 
in particular, for constructing Example l7.6.2l and for providing us with a complete 
classification of rank 2 bundles over the torus, which we hope will be useful for 
future research. 

We start from a non-trivial bundle E — > T' whose fibers are R'^"'. Such examples 
are well-known, but we detail a special one in Example 17.6.21 

As it is well-known, when 1 = 1, the only obstruction to triviality is the orien- 
tation but when / > 2, there are other obstructions to triviality. We just mention 
the Euler characteristic or characteristic classes (Whitney-Stiefel or Pontryagin for 
real bundles or Chern classes for complex bundles) . See |Ste5H IMS741 IHus94| . The 
following construction is very similar to constructions in |GS99| Section 1.4]. 

We now consider a manifold as a bundle given by 

(123) 

= r ©£:"©(r' XT'), 

where f = E, E"^ = E* — the notation E* indicates the dual bundle of linear 
functions on the fibers — and © is the Whitney sum of bundles. We use the index 
s,u to give an indication of future constructions. 

We will also introduce the notation TT' = f ^ so that we can write 

(124) 7W = © 5" © 

We denote the projections associated to each of the bundles £^,£'^,8'^ by 11*, 11", 
respectively. 

The manifold is a bundle over T' whose fibers are M'*^' x R*^^' x R'. We can 
denote points in M as (e**, e", e'^, 0), where e*^ £ {II'^)~^{9), a = s,u,c. 

We also recall that if i? is a linear bundle over a manifold JV, TE can be canon- 
ically identified as a bundle over TN with fibers isomorphic to those of E. The 
basic idea is that the tangent directions along the fibers of E can be identified with 
elements of the fibers since the space is linear. 

Hence, we will write points in T(^f,',e^,e'^.e)-M as (w*, u", w'^, i;*) where v'^ G £g , 
a = s,u,c and v* £ TeT' . Of course, we have the fact that the tangent bundle over 
the torus is trivial. 

In a coordinate patch which trivializes the bundle, we can introduce the form 
^su _ ^^"^'efdef. The key observation is that, even if the definition is in a 
coordinate patch, a change of coordinates in the patch leaves the form invariant. 
This is completely analogous to the coordinate construction of the canonical form 
in a cotangent bundle |Arn89| IGSTT) . 

We also construct the canonical one-form in £'^ by a"^ = Si=i ^idOi and consider 
the form a = a*" + a^. 

The form fl = da = d{a'''^ + a"^) is symplectic on A4. Indeed, it is clearly 
closed by definition. The fact that it is non-degenerate can be seen directly since. 
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in the coordinate patch which triviahzes the bundle, it has the standard form. As 
a consequence, M can be considered as an exact symplectic manifold. 

We now relate the previous construction to our problem. We consider a linear 
bundle isomorphism on f over a rotation T^, i.e. a family of invertible linear maps 
Ag : £§ ^e+w • We can then form a bundle isomorphism on f © over the same 
rotation which preserves the form a^" by setting 

A^(e^e") = (Aee^(A,-l)Te-). 

Then, the mapping 

F(e^e^e^^)-(A^(e^e"),e^0 + c.) 

is exact symplectic. The embedding K : T'' ^ Ai given by K{9) — (0,0, 0,0) 
clearly satisfies ([6]). If we compute the non-degeneracy conditions for this trivial 
solution, we obtain that A{6) — Id and Q{0) — Id, which is the derivative of the 
frequency on the center direction. 

The hyperbolicity condition is verified if 

P1I <Mi<l 

and 

ll(^")-'ll = 11(^^)^11 <M2<1. 

This can be arranged by multiplying by a constant if necessary. Note that in 
this case, we can take ^3 to be as close to 1 as desired. 

Furthermore, if G is analytically close to F (i.e. \\F — G\\b < £, where S is a 
suitable complex subset of A4) and exact symplectic, then we have 

\\GoK-KoTjp, = \\FoK -GoK\\p„ <e 

so that if £ is small enough the hypotheses of Theorem 13.111 are met. 

7.6.2. An explicit example. To make the whole construction more concrete, we just 
end with an explicit example of a non-trivial R^-bundle over T'^ with positive Euler 
characteristic explained to us by Prof. Gompf. Many more examples can be found 
in |MS74| . Applying the construction in this section to these examples gives us sym- 
plectic manifolds and whiskered tori with non-trivial stable and unstable bundles. 
We construct a bundle over §2 with non-zero Eulcr characteristic. If we identify 
with C using the standard identification and §^ with the Riemann sphere, we 
can construct a non-trivial bundle in the semi-sphere, whose boundary is the circle 
S"'^ = {\z\ = 1}, by identifying the product bundle. We just give a gluing map on 
the unit sphere bundle, and extend it homogeneously. Hence, it suffices to give an 
identification mapping i from §^ x to itself. The first factor is the boundary 
of the disk and the other factor is the unit bundle. We take i(z, w) = (1/z, z"w). 
Using partitions of identity, one can extend this bundle on a disk to a bundle of 
the torus. 

8. Finite-dimensional flows 

This section is devoted to the application of our method to find invariant tori for 
symplectic (locally Hamiltonian) vector-fields. Although we have already presented 
a result — in a rather abstract way — on existence of invariant tori for vector-fields 
in Theorem 13.151 we now present a direct proof of the results, following similar 
methods as in the case for maps. One motivation for writing this section is that 
the proof leads immediately to algorithms, which may be useful for applications 
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involving vector-fields rather than maps. It may be of interest for practitioners to 
have algorithms for flows. 

The proof for flows can also serve as a starting point for a proof for PDE's. 
We also note that the methods developed here apply to some ill-posed partial 
differential equations, which do not admit time-1 maps. Of course, the adaptation 
of the strategy of proof to PDE's involves several technical considerations (the 
generators of the evolutions are unbounded operators rather than differentiable 
ones). We postpone these considerations on PDEs to a forthcoming paper (see 
[dlLSOTp . 

We will study first the case of locally Hamiltonian flows. The case of globally 
Hamiltonian flows will be discussed in Section [9l 

8.1. Some preliminaries on symplectic geometry. In this section we recall 
several well-known facts on symplectic geometry of vector-fields. 

We will consider vector-fields on an exact symplectic manifold A4 with symplectic 
structure = da. We have the following definitions. 

Definition 8.1. We say that a vector-field X on M is symplectic when 
where Cx stands for the Lie derivative with respect to X . 

Definition 8.2. We say that a vector-field X is exact symplectic when there exists 
a smooth function W on M such that 

Cxa^ dW. 

An easy calculation checks that exact symplectic vector-fields are symplectic: 

Cx^ = Cxda = d{Cxa) = d{dW) = 0. 

However, the converse is not true. A well-known example is the following: consider 
the manifold = T x R. We denote the corresponding coordinates {q,p) and we 
set a = pdq and Q — dpAdq. Consider now the vector-field X — dp. It is symplectic 
but not exact symplectic. 

Using Cartan's formula and the fact that dil = 0, we obtain that X is symplectic 
if and only if 

(125) dix^ + ixdVL = dix^- 

This means, by Poincare lemma, that locally we can write 

ix^ = dH. 

Of course, (|125p does not imply that H is a. global function since in general it is 
only locally defined. 

As a matter of fact, H will be a global function if and only if the vector-field X 
is exact symplectic. Indeed, since 

dW — Cxct — d{ixce) + ix^, 

we see that, if X is exact symplectic, we can take H = W — ixct. 

The above discussion shows that the only difference between symplectic and 
exact symplectic is the (de Rham) cohomology class of ix^. We introduce the 
following definition. 
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Definition 8.3. Let K be an embedding from T' into A4. We say that a family of 
vector-fields X\ with A G M' spans the cohomology of K{T^) at X ^ X if the map 

^ ^ lx[^*^x,^l]ix^jv 

is an isomorphism. Here we denote i?^(T') the first de Rham cohomology group of 
T', which is well known to be M' (see |Hat02] ). 

In T' X R' with the standard symplectic form, we have that, denoting by pi the 
coordinates along E' , the family 

; 

i=l 

spans the cohomology at every A. Of course, in this case, the cohomology classes 
have a very simple characterization as the averages along each of the elementary 
cycles of T' . 

8.2. Setting of the equations. The result for flows is based on the study of the 
equation 

(126) d^K{e) = x{K{e)), 

for K : Dp D T' M, where the operator d^^ (derivative in the direction oj) is 
defined by 

i—l 

and the vector- field X : A4 ^ TAi is symplectic and real analytic. 
Let St be the flow of X. If if : T' ^ is a solution of (|T2H1) then 

(127) St{K{e)) = K{e + Lut), 61 gT', teR, 

and therefore the range of K is invariant by St- Indeed, considering 6' e T' fixed, 
both sides of (|127p satisfy the same Cauchy problem. 

We first deal with a family of vector-fields X\ and we prove a version of the 
translated torus theorem. For an exact symplectic vector-field we will embed it 
into a family, then prove a vanishing lemma and finally prove the existence of an 
invariant torus. For families X\, the equation under consideration is 

(128) d^K{e) = Xx{K(e)), 

where A G M', the dependence of X\ in A is at least and we assume that the 
vector-field X\ spans the cohomology of Ko{T'') in the sense of Definition l8.3i where 
Kq is an approximate solution of (|128p . 

A very important role will be played by the linearized equation 

(129) ^=Ax{e + cjt)A, 

where A\{6) = DX\{K{6)). Since Ax is a bounded operator, equation (|f 29p admits 
an evolution operator, which is defined for all < G M, and we will denote it Ug{t). 
It is characterized by 

(130) jUg{t) =Ax{e + ujt)Ue{t), Ue{0) = Id. 
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8.3. Non-degeneracy conditions. To establish the existence of tori, we will re- 
quire non-degeneracy conditions similar to the ones considered in the case of maps: 
namely, a spectral condition and a twist condition. 

Condition 8.4. (Spectral non- degeneracy condition) Given X gM} and an embed- 
ding K : Dp D ^ M we say that the pair (A, K) is hyperbolic non- degenerate 
for the functional equation p28p if there is an analytic splitting 

TK(e)M = Sxie) ® ^K{e) ® ^K(e) 
invariant under the linearized equation (I129P in the sense that 

Moreover the center subspace f^^^-j has dimension 21. We denote 11^^^^, ^Kie) '^^^ 
projections associated to this splitting and we denote 

U^'''-(t)^Ue(t)\,^^^. 

Furthermore, we assume that there exist /3i , /32 , /Ss > and Cu > independent 
of 9 satisfying f]^ < Pi, (3^ < P2 o,nd such that the splitting is characterized by the 
following rate conditions: 

\\UI{t)UI{T)-X < Cue-P-'^'-^\ t>r>0, 
(131) \\Unt)U^{r)-'\\p < CueP^-^'-^\ t < r < 0, 

\mm^ir)-'\\p < ae^^l*-^!, t,reR. 

Remark 8.5. As in the case if maps, if we have an approximately invariant splitting 
and 

||C/|(t)f/|(r)-i||, < e-^^^'-^\ T/2<t-T<T, 

\\U^m^{r)-'\\p < e^^(*-^), T/2<T~t<T, 

\mm^e{r)-'\\p < e^^'l*-^! T/2<\t-T\<T, 

for some T large enough, then there exists a true invariant splitting, close to the 
approximately invariant one, and the bounds (|13ip with respect to this new splitting 
hold. This can be checked by using the time T map. 

Remark 8.6. The previous non-degeneracy condition just expresses that we can as- 
sociate semi-groups in positive and negative times to the operator Ax(9-\-ujt). More 
precisely, since the systems under consideration are non-autonomous, we should 
write 

f = ^, m = 9. 

Note that if the systems were autonomous, the exponential bounds would follow 
from the spectral properties of A\. 

The linear operators Ug''^'^{t) enjoy the following co-cycle property. 

Lemma 8.7. For all 9 and lu and all times t,T we have 
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Proof. It follows from the classical argument of uniqueness for Cauchy ODE prob- 
lems. Dropping the indexes s, c and u, for 9, lo and t fixed, we define the functions 

^i,t(T) = Ueit + r)7/-o, V2,t(T) = C/e+^t(T)C/0(i)Vo 

for an arbitrary -0o- Since J76»(0) is the identity operator, these two functions satisfy 
the same Cauchy problem and hence are equal. □ 

Condition 8.8. (Twist condition) Let Ax(9) = DXx{K{6)) and 

N{9) = [DK{9)'^DK{e)]-\ 

We say that the pair (A, K) satisfies the twist condition if the average on T' of the 
matrix 

Sx{9) = N{9)DK{e)'^[d^{J{K)-^DKN) - AxJ{K)-^DK N]{9). 
is non-singular. 

If a pair {\,K) with K : DpZ)T^ ^ M satisfy both Conditions [Ol and [8?8l we 
write (A, K) e ND{p). If X does not depend on A we simply write K e ND{p). 

We note that Conditions 18.41 and 18.81 hold in open sets of K. The fact that 
Condition 18.81 holds for an open set (in the topology) is obvious since it is 
the non-degeracy of a matrix that is just an explicit algebraic expression involving 
derivatives. The fact that Condition 18.41 is stable under perturbations will be the 
content of Section 18.71 

8.4. Statement of the results. The first result below provides an existence re- 
sult in the case of a family of symplectic vector-fields. From a sufficiently good 
approximate torus for a vector-field in the family it provides an invariant torus for 
a translated (with respect to the parameter) vector-field in the family. 

Theorem 8.9. Let lu G DhiKjv) for some k > and v > I ~ 1. Assume the 
following hypotheses 

(1) The vector-fields Xx are symplectic for every A G M'. 

(2) The family Xx spans the cohomology of Kq{T^) at X ^ Xq in the sense of 
Definition \8.3[ 

(3) The pair [Xa^K^) satisfies the non- degeneracy Conditions \8.4\ and \8.8i 

(4) The vector-fields Xx are real analytic and they can be extended holomor- 
phically to a complex neighborhood of the image under Kg of Dp^ : 

Br = {ze <C^'^\ 39 e {|Im6l| < po} s.t. \z - Kq{9)\ < r} , 

for some r > 0, and are with respect to X. 
Define the error Eq by 

Eoi9) ^ dMS) - Xx,iKoi9)). 
Then there exists a constant C > depending on I, v, \Xx\c'^(Br)j ll-C^o||poj 
PO' II OA ''llpoj ll'5'ollpo! l(8'Vg(<5'o)) ^1, (where So and Nq are as in Condition 
\8.^ replacing X by Aq and K by Kq) and the norms of the projections \\^%'^g(g-f\\po 
such that, if Eq satisfies the estimates 

Ck^S-^''\\Eo\\p„ < 1 

and 

C^5-"'\\EoU<r, 
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where < (5 < min(l,po/12) is fixed, there exists an embedding Koo and a vector 
Aoo G such that (Aoo, Koo) G ND{poc '■= po — 6(5) and 

(132) d^KooiO) ^ Xx^iKooie)). 
Furthermore, we have the estimates 

ll^oo-ifollp^ <Ck^S-^''\\Eo\\p, 

and 

\Xoo-Xo\<Ck^S-^''\\Eo\\p„. 

The following theorem deals with the existence of invariant tori for exact syni- 
plectic vector-fields. It follows from the translated torus version Theorem 18.91 ap- 
plied to a suitably chosen perturbation of the exact symplectic vector-field X and 
a vanishing theorem whose proof is postponed to Section 18.81 

Theorem 8.10. Let lu e Dh{n, i') for some k > and v >l — 1. Assume that 

(1) The vector-field X is exact symplectic. 

(2) Kq satisfies the non-degeneracy Conditions \8.4\ and \8.8l 

(3) The vector-field X is real analytic and it can be extended holomorphically 
to a complex neighborhood of the image under Kq of Dp^ : 

Br ^{ze C^'^l 36 e {\Im9\ < po} s.t. \z - Ko{0)\ < r} , 

for some r > 0. 

Denoting Eq the initial error, there exists a constant C > depending on I, u, 
\X\c'^(B,.), \\DKo\\p^, \\No\\p„, llS'ollpo, l(avg(S'o))"^|, (where Sq and Nq are as in 
Condition l 8. 8\ revlacina K by Kq) and the norms of the projections \\^'Kg{e)\\po such 
that, if Eq satisfies the estimates 

Ck^S-^'^WEqWp,, < 1 

and 

CK^S-^-'WEQWp^Kr, 

where < S < min(l,po/12) is fixed, then there exists an embedding Koo G 
ND{poo '■— Po — 6^) such that 

(133) d^Koo{e)=X{Koom- 
Furthermore, we have the estimate 

\\Koo-Kq\\p^ <CkH-^''\\Eq\\p,. 

Remark 8.11. One could also formulate a local uniqueness result in the case of 
vector-fields. This can be done by a reduction to a time-one map (see |Dou82j l. 

8.5. Linearized equation. In this context we define the operator 

g^{X,K)^d^K -Xy^oK 

and we want to solve the equation Qu){X,K) = 0. As in the case of maps this 
will be done through a KAM iterative procedure, starting with (Aq, Kq) such that 
E ^ Guj (Ao, i^o) is sufficiently small. Therefore we are lead to consider the linearized 
equation 

(134) d^A{e) ~ A^i9)A{e) - ^-^A^^W) A = -E{e), 
where Ax{e) = DXx{K{9)). 
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Let ^ : T' ^ be a function. From the spectral non-degeneracy condition we 
have 

(135) TiKie+^t)Ue{t)m = UemK^g^m, 

where 11 stands for any of the projections H'^ and H". Differentiating with 
respect to t both sides of (|135p and using (|130p we obtain 

= A^{9 + u;t)Ug{t)UKie)m- 
Evaluating this expression at i = and using the definition of di^ we have 

which implies 

(136) UKie) [5^ - A^idMO) = [9. - Ax{d)]UKie)m- 

8.5.1. Linearized equation on the center subspace. We first project equation (|134|) 
on the center subspace. Using (|136p we immediately obtain 

(137) d^A%9) - Ax{e)A^9) - nj^^.^^M^A = -E%9), 
where A%e) = n^(g)A(6') and £"^{9) = n^(g)£;(0). 

8.5.2. Small divisors equations and isotropic character of the torus. The following 
result, which is completely analogous to Proposition 14.41 deals with the resolution 
of small divisors equations along characteristics (see |Riis76a| . |Riis76b| . [Riis75| . 
(dlLOlcQ . 



Proposition 8.12. Assume that oj £ Dh{K,v) with k > and v > I — \. Let 



h : Dp D T'' ^ Ai be a real analytic function with zero average. Then, for any 



< S < p there exists a unique analytic solution v : Dp^s Z) Ai of the linear 

equation 

j=i J 

having zero average. Moreover, if h E Ap then v satisfies the following estimate 

\\v\\p^s<Cn5-''\\h\\p, Q<5<p. 
The constant C depends on v and the dimension of the torus I. 

The following result provides the approximate isotropic character of the torus. 
This proposition is similar to the one in |dlLGJV65] and we do not reproduce its 
proof here. We note that it also follows by taking time-1 maps from the corre- 
sponding result for maps, which we have established in Section B. 1.11 

Proposition 8.13. Let K : Dp D T' — > A^, p > 0, be a real analytic mapping. 
Define the error 

E{9) ■.= d^K{9)-Xx{K{9)). 

Let L{9) = DK{9)^ J{K{9))DK{9). There exists a constant C depending on I, v 
and \DK\p such that 

\\L\\p-2i<CK8-^-^^\\E\\p, 0<5<p/2. 
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Once again, we use a normalization argument which aUows us to write equation 
()137|) in a suitable form. To do so, we need a result which allows to approximate 
the center subspace with the range of the 2d x 2Z-matrix 

(138) M{9) = [DK{e), J{K{9))-^DK{e)N{9)], 

where N{9) is the normalization / x /-matrix given by N{9) = [DK{d)^ DK{9)]~^ , 
as in Proposition 14.161 One can prove the following result. 

Proposition 8.14. Denote hyV}i(^g-^ the range of M {9) andbyH^^g^ the projection 
onto ^K{e) according to the splitting S^f^g-^ © ^k(0) © ^k{0) ■ 
Then there exists a constant C > such that if 

S-'\\E\l<C 

then we have the estimates (here distp has to be understood as the distance of 
subspaces in the Grassmanian sense) 

distp_25(r^(e),f^(e)) < C5-^\\E\\p, 

(139) „ , 
m'k(e)-^lie)\\p-25<C6-^\\E\\p 

for every 5 £ (0,p/2) and where C, as usual, depends on the non-degeneracy con- 
stants of the problem. 

The proof of the previous proposition follows the same lines as the one of Propo- 
sition HTTBl We refer the reader to CoroUarv lS. 221 where we construct exact invariant 
splittings from approximate ones. 

We introduce the change of function A'^ = + e^, where ^ : T' ^ TAA, with 
C(^) e TKie)M and e = n^^g^ - U^j^^^y We then get 

(140) [djd{9) - A,,{9)M{9)]m + M{9)d^m - n^(,)^^^||^A = -E^{9), 

where we have dropped the terms depending on e^, which are quadratic in the 
error. As in the case of maps, the matrix M{9) is not invertible but the matrix 
M{9y J{K{9))M{9) is. Multiplying equation (fTiPl) by M{9y J{K{e)) and then 
by [M^ J{K)M)~''^ , we get the following equation 

{M{ey J{K{9))M{9))-Hl{9y J{K{9))[dJl{9) - A^{e)M{9)]m + 9^9) 

= (M(g)^J(g(g))Af(g))-^M(g)^J(g(g))[n-^(,) ^^^^^^^^U - E^{9)]. 

We are going to normalize the matrix d^M{9) — Ax{9)M{9). To avoid some compu- 
tational technicalities, we perform this normalization only when K is a, solution of 
()128p . We refer the reader to the case of maps on how to handle the computations 
in the approximate case. 

Lemma 8.15. Let {X,K) he a solution of 

(141) d^K{9) = Xx{K{9)) 

and M be the matrix defined by (jl38p . Then there exists a I x l-matrix Sx{9) such 
that 

(142) d^M{9) - Axi9)M{9) = Mi9) (^^'^ ^^P) ' 
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The matrix S\{0) has the form 

Sx{e) = N{9)DK{ey[d^{J{K)-^DKN) - AxJ{K)-^DK N]{e). 

Proof. Exactly in the same way as in the case of maps, if is a solution of (|141|) 
the columns of M generate the center subspace. Since — A\{9) commute with 
^'k(0) have that 

(143) d^M{e) - Ax{9)M{9) = M{e)C{e) 

for some 21 x 2/ matrix C{6). Differentiating equation (|14ip with respect to 9 we 
obtain 

(144) d^DK{e) = Ax{e)DK{e). 
This implies that 




Identifying blocks in (|143p we end up with 

(145) d^{J{Ky^DK N) - AxJ{K)-^DK N = DKSx + J{K)~^DK N Rx. 

Multiplying (|145p by DK^ J{K) and using the isotropic character of the invariant 
torus, i.e. 

L{e) = DK{ey j(K{e))DK{e) = o, 

it follows that 

(146) Rx =DK^ J{K)[d^{J{K)-^DK N) ~ AxJ{Ky^DK N\. 
Expanding d^{J{K)^^ DK N) and using equation (|144p . we get 

d^{J{K)-^DKN) = d^{J{K)-^)DK N + J{K)-^AxDK N + J{K)-^DKd^N. 

By differentiation of NN^^ ^ Id, using ()144l) we easily obtain 

d^N = -NDK^[Aj + Ax]DKN. 

Also d^{J{K)-^) = -J{Ky^DJ{K)AxDK J{K)-\ 

Moreover the symplectic character of the vector-fields Xx, i.e. Cx^,^ = can be 
expressed by (recalling the definition of the Lie derivative) 

(147) ^[D$7j(<i>0i?$t]|t=o = 0, 

where $4 is the flow solution of Xx and (|147p implies 

Ajj{K) + J{K)Ax + DJ{K)X{K) = 0. 

Using the previous calculations we obtain that the right-hand side of (|146|) vanishes, 
i.e. Rx = 0. 

Now multiplying (|145p by N DK^ and using the definition of N we have 

(148) Sx{d) = N{e)DK{ey[d^{J{K)-^DKN) - AxJ{K)-^ DK N]{9). 
Using again the previous calculations we can express Sx as 

Sx = N DK^ J{Ky^\id2d - DKN DK^]{Ax+ Al)DKN. 

We emphasize that this last formula coincides with (I148P only when K is an exact 
solution. If K is only an approximate solution then both expressions are approxi- 
mately equal. 

□ 
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We now turn to the case of approximate solutions. The procedure is similar to 
the one of the case of maps. 

When K is just an approximate solution, we define 

(ei, 62) = d^MiO) - A^{0)M{9) - M{0) (^[J^ ^"^^^^ . 

Some computations, using that d^DK{9) — A\{d)DK{9) = E{9) and the dcfintion 
of Sx give ei = DE and 62 = 0{\\E\\p, \\DE\\p). 

Next we just state the result without proof, but we indicate that it is quite 
analogous to the proof in the map case. We first identify - up to a small error - 
the center space with the span of the tangent and its symplectically conjugate and 
then compute the matrix of the derivative in these coordinates. 

Lemma 8.16. Assume uj e Dh{n, i') with k > and v > I — \ and \\E\\p is small 
enough. Then there exist a matrix B{9) and vectors pi and p2 such that equation 
(|140p can be written as 

[jj +B{9)]m+dM0)^Pi{0)+P2{e) 

(149) -(M(0)Tj(^(0))M(0))-iM(0)Tj(i^(0))n^(,)^^^^|MA. 

Moreover, the following estimates hold: 

(150) lbi||p<C||i?||p, 

where C just depends on \\J(K)\\p, ||A^||p, H-Dii'llp and ||n^^gj||p. For p2 and B we 
have 

(151) \\p2\1-2S<Ck5-^'^+'^\\E\\1 
and 

(152) \\B\\p_2S<CkS-^''+^^\\E\\p 

for 5 G (0,p/2), where C depends on I, v, ||A^||p, ||-Diir||p, \Xx\(j2(^b^^, \ J\c^[B.r) '^'^^ 
l|llA'(e)llp- 

8.5.3. Solution of the reduced equations. The solution of the reduced equations 
works in the same way as in the case of maps. We sketch the procedure in this 
section and we emphasize on the cohomology obstructions on the equations. 
We write ^ = (^1,^2). We introduce the operator 



(153) CS, 



(0! ^Q^)^^d^^^PM + Qm. 



where pi — (pii,pi2) and Q = {Qi, Q2)- Using this decomposition of f^^g-j we can 
write equation (jl53p in the form 

S{e)UO) + dM{e) = pn{9) + Qi{9)A, 
d^US) = Pi2{0) + Q2{0)A- 

We furthermore have 

= (iV^iJA-^J(X)-^)(g)[(jJXiVJj^^)(g)-Id2rf]J(i^(g))n^(g) ^^"™^ A, 



Q2{9) = i?if(0)Tj(if(0))n^(,)^M|MA. 
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The assumption of spanning the cohomology of K{T'') for the vector-field X\ en- 
sures that we can choose A such that tlic second equation is solvable in the sense 
of Proposition 18. 121 Indeed, notice first that the cohomology along the hyperbolic 
bundle of the form K*ix^^ is trivial and we have then 

Identifying the cohomology class of a form in i?^(T') to its integral on the torus T' 
and using the fact that the family X\ spans the cohomology of K{T^) at A gives the 
result (since we can choose A such that the average of pi2{0) + Q2(^)A vanishes). 

The degree of freedom we get on the average of ^2 then allows us to solve the 
equation on ^1. Recall that we use the non-degeneracy conditions as stated in 
Condition 18.81 We obtain the following proposition. 



Proposition 8.17. Assume ut E Dh{K, v) with k > and v >l ~1, and {X, K) is 
a non-degenerate pair. If the error \\E\\p is small enough, there exists a mapping 
^, analytic on Dp_25 and a vector A G M' solving equation (|153p . 

Moreover there exists a constant C > depending on v,l, \\K\\p, |(avg (A))~^|, 
||A^||p and W^'j^^g^Wp such that 

and 

\K\<C\\E\\p. 

8.6. Linearized equation on the hyperbolic space. We project the linearized 
equation (|134p 

on the stable and unstable subspaces by using the projections H^j-g^ and H^j-g^ 

respectively. We denote A'* (6*) = n^(g)A(6l), A" (6*) = n^(g)A(6l) and E{e,\K) = 

^^^4piA -£((?). 

Using the previous notation and (|136p we obtain 

(154) d^^'{e) - Ax{e)/\'^{e) = n'K(e)E{e. a, a) 

for the stable part and 

(155) a^A"(0) - Ax{e)^^{d) = n^(,)^(0, a, a) 

for the unstable one. 

The following result provides the solution of the previous equations. 



Proposition 8.18. Given p > 0, equations ()154p and ()155p admit unique analytic 
solutions A'* : Dp f and A" : Dp £^ respectively, such that A*'"((?) G ^i^gy 
Furthermore there exist constants C"'" such that 

(156) ||A^^"||p<C^^"(||i?||p + |A|), 

where C^"^ depend on (3i, ||n^(g)||p (resp. P2, \\^Kie)\\p) and Ch, \\ ^^g^^^ ||p- 

Proof. The proof is based on the integration of the equation along the characteris- 
tics 6 + ujt and the use of the spectral non-degeneracy Condition 18.41 We give the 
proof for the stable case, the unstable case being symmetric (for negative times). 
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We introduce the function A{t) = A^iO + ujt). If A"* has to satisfy (fTSi]) then 
A{t) has to satisfy the equation 

(157) ^ A(t) - Aa(0 + u;t)K{t) = n^Kie+^,)E{e + u;t, A, A). 

We first derive heuristicaUy a formula p60p for A'*. Then, examining the for- 
mula, it will be easy to justify the derivation. 

Let Ug{t) be the evolution operator characterized by 

(158) ^Ue{t) - A)^{9 + uot)Ue{t), Ue{0) = Id. 
at 

Using the formula of the variation of parameters we have 

(159) A{t) ^ Ueit) [a(0) + /* U,\s)Wj,^,^^^)Eid + us, A, A) ds 

Jo 

Using the co-cycle property given by Lemma [8771 we have Ug^{s) = C/e+[js(~s)- 

Since formula (|159p is valid for all 9 G Dp D T' we can use it substituting 9 by 
9 — Lot and recovering the notation A*': 



L Jo 

We assume that A*, the solution we are looking for, stays in and it is bounded; 
then Ug^^t{t)A^{9 — cut) goes to when t goes to oo. Using again the co-cycle 
property we have 

Ug-ut{t)Ue-^{t~s){-s) = Ue-u>{t-s){t ~ s). 

Then we write 



A' {9) = Ue-Ut)^''{0~^t)+ / Ue-u,it-s){t-s)n'j,^,,^_..E{9-u:{t~s),\,K)ds. 

Jo 

Performing the change of variable t — t — s and letting i go to cxd wc finally obtain 

/>OC 

(160) A^(0)=/ Ue-ur{rWK(9-u.r)E{0-LOT,\K)dT. 

Jo 

Using the spectral non-degeneracy hypothesis the subintegral function is bounded 
by C/ie~''i'^||n^^g_^^^ llplji^llp. The exponential bound assures the convergence of 
the integral and also permits to obtain the bound 

l|A'^l|p = (c,//3i)||n^(,__)||p||^||p. 

Once we have formula (|160p we check directly that A* is indeed a solution of (|154p . 
The absolute convergence of (|160l) justifies the exchange of limits and rearrange- 
ments used in the derivation. 

The uniqueness follows from the fact that if we start with any solution Af of 
P54p . doing the previous manipulations we will end up with the same explicit 
formula p60)) . □ 
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8.7. Change of non- degeneracy conditions in the iterative step. The next 
result deals with the measure of the change of the splitting when we perturb a 
linear system in an Euclidean space M. 

Let Ax{0) be a family of linear maps from an Euclidean space M into iteself, 
depending on 6 G Dp D T'^ and A G M' and let 11$ be its evolution operator, i.e. 

j^Ueit) = Ax{e + cot)Ue{t), Ue{0) = Id. 

Assume that M has an analytic family of splittings 

M = £^® £^ © £^ 

invariant by Ue in the sense that Ue{t)£g''^'^ = £g'^'^t- H^''^'^ the projections 
associated to this splitting and Ug''^'^{t) = = Assume furthermore there 

exist (ii, P2, > and C/j > independent of 9 satisfying P3 < Pi, < P2 and 
such that the splitting is characterized by the following rate conditions: 

mm^irr'Wp < Che-^'^'-^\ t>T>o, 

WU^m^ry'Wp < Che^^^'-^\ t<T<0, 

\\mMir)-'\\p < Che^'^'-^\ i,T e K. 

Proposition 8.19. Assume that ^a(^) *s a family of linear maps as before. Let 
A\{6) he another family such that \\A\ — A\\\p is small enough. Let Uo{t) denote 
the evolution operator corresponding to Ax, i.e. 

jUe{t)^Ax{e + u:t)Ue{t), [/^(O) = Id. 
Then there exists a family of analytic splittings 

M = £l® £e ® £0 
which is invariant under the linearized equation 

^A = Ax{9 + u;t)A 
at 

in the sense that 

We denote n^'*^'" the projections associated to this splitting and denote 

Then there exist Pi, P2, Ps > and Ch > independent of 6 satisfying Ps < Pi, 
Pa < P2 o-nd such that the splitting is characterized by the following rate conditions: 

mntWeirr'Wp < Che^'^'-^\ t<T<0, 

Wmtmirr'Wp < Che^'^'--\, t,reR. 
Furthermore the following estimates hold 

(161) \\Ill'''^-Ul''''^p<C\\Ax-Ax\\p, 

(162) \pi-pi\<C\\Ax-Ax\\p, i = 1,2,3, 

(163) Ch = Ch. 
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Proof. We provide the proof of the statements concernmg the stable subspace. We 
divide it into several steps. We use the notation of Condition 18.41 

Step 1: Construction of the invariant splitting. We look for the invariant 
splitting associated to the linearized equation 

(164) ^W(t) ^ Ax(e + ujt)W(t) 

at 

focusing on the stable bundle. We write (|164|) as 

(165) ^^"^^^ = AxiO + Lut)W{t) + Bx{9 + ujt)W{t) 

with Bx = Ax — Ax- Since we are interested in solutions decreasing exponentially 
at oo, for a > we introduce the space 

Ca^ {f ■■ [0, oo) ^ C^^^ I / continuous, supe"*|/(OI < oo}, 

t>o 

with the norm \f\a = supt>o e°*|/(i)|. 

Given a G (/33, /9i) we look for solutions of (jl65p in the space Cq. We begin with 
the following auxiliary result. 

Lemma 8.20. Let 6 G Dp D T', ^ G £^(^g-^ and H e Cq, with, a G (/33,/?i). Consider 
the equation 

(166) w' ^ Ax{e + ujt)w + H{t). 

Then there exists a unique Junction 1C(S^,H) G Cq, such that 

(i) IC{^,H) is solution of ()166p . 

(ii) n|/c(c,7?)(o) = e 

Moreover IC{£,,H) = /Ci(^) + IC2{H), where /Ci : ^1 — > Cq and /C2 : Cq — > Cq are 
bounded linear operators and 

(167) \ICi\<Ch, 

(168) |/C2|<cJ-^ + ^i 
where = sup,g^jn^'™|. 

Proof. If w G Cq is a solution of (|166p in [0, 00) and t, r > we have 

(169) w{t) = Ue{t)Ue{T)-^w{T) + Ug{t)Ue{s)-^H{s) ds. 

Projecting (|169p to the center-unstable subspace and using the invariance of the 
splitting £g {£g £g ) with respect to Ug and writing 11™ the projection onto 

CC fT. cu 

On ki) On 



(170) n™^,«;(i) = c/e(i)c/e(T)-^n™^,w(T) + j Ug{t)Ug{sr'nT+^Ms)ds. 

If r > i we have 

\Ue{t)Ue{T)-'llT+.Mr)\ < C,e^^'(-"*)|n§!;_| e-^^^^U 
which goes to zero as r tends to 00. Also, if s > t 

\Ug{t)Ug{s)-'lir+^,H{s)\ < C,,e^^^^-'^\nr+^Je-^\H\^ 
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guarantees that we can take limit t oo hi the integral in ()170|) . Then we have 

Kl.Mt) = f Ue{t)Ue{s)-'llT+^Ms)ds. 
Using the projection to the stable subspace, we obtain 

w{t) =ni^^Mt) + nr+^Mt) 

(171) =C/e(i)nXO)+ / Ueit)Ug{s)-^n'g+^^H{s)ds 

rt 







/ Ue{t)Ue{s)-^TiT+^,H{s)ds. 

J oo 



Once we have the explicit expression of w, we easily check that it actually belongs 
to Ca- We define /Ci(^)(i) = Ue{t)(, and IC2{H){t) to be the sum of the two integrals 
in ()171|) . A simple calculation gives the bounds p67p and (I168p . □ 

By Lemma [8.201 the solutions of p65|) belonging to Ca satisfy 

wit) = /Ci(n^w(o)) + )C2{Bx{e + ij-)w){t). 

Note that B\{d + uot) is bounded in t and moreover \B\{9 + u)t)\ < 7, where 
7 = II — ^aIIp- We introduce the linear map IC2 ■ Ca ~^ Ca defined by 

iC2{w)^IC2{Bx{e + uj-)w). 

Clearly \IC2{w)\ < \IC2\\Bx{0 + uj-)\. With the above introduced notation, given 
€ £g, there is a unique solution w E Ca such that ligw{0) = ^ which is given by 

w = /Ci(0 + iC2{w). 
Since |i?A| < 7 < 1 we can write 

W^{ld-t2)-^ICi{0. 

Therefore is the graph of 

00 

C ^ M^m ng"(Id - ^2)-^/Ci(e)(0) = nr J2 ^2^1 (0(0), 

fc=l 

where the sum stars with fc = 1 because II^/Ci — 0. Note that the analyticity in 9 
is preserved in all the previous manipulations, hence Mg depends analytically in 9. 
Since I/C2I < C7 then ||M'*(0)||p < C7. In a completely analogous way we find 
and integrating with negative times we get Eg and Eg''. Finally Eg = Sg'^ fl £g". 

Step 2. Estimates on the projections. To get the bounds for the projections 
we follow the same argument as in the case of maps. We only give the argument 
for the stable subspace. Let M™(0) be the linear map whose graph gives fg". 

We write 

n^e = (r,o), ui^ = {e,M%9)t), 

n™^ = (0,^™), ng^c = (M™(6')|™,|™), 
= r + M™(6')|™, 



and then 
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Since M^{9) and M'^"[9) are 0(7) we can write 
t \ _ ( Id Af™(6l) 

y M''(6') Id ; 

and then deduce that 

step 3. Estimates on the growth conditions. 

To get the exponential bounds let 

i,{t)^Ue{t)iJe{T)-^^{T) 

with iI]{t) = M{6 + ti'T)^) G Sg^^T-- The function -0 satisfies equation (|165p and 
hence 

|V(t)| < \Ue{t)Ue{Tr^^{T)\+ \Ug{t)Ug{s)-\A^ - Ax){e + ujs)ij{s)\ds, 

for t > T. Let X be the auxiliary function defined by = e'^i*|V^(<)|. Using the 
bounds of Condition 18.41 we have 

X{t) < Chxir) + ChC^ ^ x{s) ds, t > T. 
By Gronwall's lemma we have x{t) < Chx(T)e'^''^'^(*""'' and hence 

We conclude that 

|(7e(i);7e(T)-V(T)| < C,e-('5i-^'^C'7)(t-r)|^(^)|^ t > ^. 

We take Cu = Ch and /3i = /3i — ChCj, which proves (|162p . 

□ 

The first consequence of Proposition 18.191 is that in the iterative step the small 
change of K produces a small change in the invariant splitting and in the hyper- 
bolicity constants. 

Corollary 8.21. Assume that (A, K) satisfies the hyperbolic non- degeneracy Con- 
dition \K4\ and that \\K — K\\p is small enough. If we denote A\(6) = DX\{K{6)), 
we can define an evolution operator, denoted Ugit) such that 

^Ue{t) = A^{e + u:t)Ue{t), f/e(0) = Id. 
dt 

Then there exists an analytic splitting for K , i. e. 

which is invariant under the linearized equation (jl29p (replacing K by K) in the 
sense that 

We denote 11^^^^, ^'k(e) ^^'^ ^^k[e) projections associated to this splitting. De- 
noting 

K {8) 



68 



E. FONTICH, R. DE LA LEAVE, AND Y. SIRE 



there exist /?i , /32 , /?3 > and Ch > independent of 9 satisfying < f3i, < (32 
and such that the splitting is characterized by the following rate conditions: 

WUimi^rr^Wp < C,e-^^(*-^), < > T > 0, 

\\Unt)U^{T)-^\\, < C,e^^(*--), < < T < 0, 

\\mtmir)-'\\p<Che^'^'-^\ t,reM. 
Furthermore the following estimates hold 

(172) \\U%^---U^^^^]\\,<C\\K~K\\,, 

(173) \/3,-P,\<C\\K~K\\p, * = 1,2,3, 

(174) Ch = Ch. 

Proof We just take Ax{0) = DXx{K{9)), Ax{e) = DX^{K{9)), f^^^"; = 

R{e) ~ ' ^^K{e) ^ ^^e ^(^) ~ ^ "^^ Proposition 18. 191 and we use 

that \\Ax{6) - Ax{e)\p < \\X\\c2\\K{e) - K{e)\\p. a 

The second consequence of Proposition [STTH is that if we have a sufficiently good 
approximate splitting associated to equation (jl30p then there is a true invariant 
splitting nearby. 

Corollary 8.22. Assume that Tk(9)M — ^^^(e) ® ^K(e) ® ^'^ splitting ap- 

proximately invariant under the linearized equation (|129p with evolution operator 
Ue{t), in the sense that A\(9) — DX\{K{9)) can he represented as 

(A\\9) Af{9) Af{9)\ 
A,(0)= \Af{9) Af{9) Af{9)\ 
\Af{9) Af{9) Af{9)l 

with respect to this splitting with ||^y(^)||p < CJ^"'" ||i?||p if i ^ j. Letli*^'^^^ he the 
projections associated to to this splitting. 

Let I7g'°'" he the evolution operators of A = A\^{9 + ujt)A, A = Af{9 + ujt)A 
and A = A^^19 + tot) A respectively, and assume 

\\UI(t)UI{T)-X < Cle-^'^^'-^\ t>T>0, 
\mt)U^{r)-'\\p < Cle^'^^'-^\ t<r<0^ 
\\mmr)-^,<C:e^'^\'--\, t,TeR, 

for some /?i 2 31 C*;* > such that < PI, P^ < Pt- Then there exists an analytic 
splitting Tx(e)M^ — ^K{e)®^K(e)®^K(e) "invariant under equation (|129p . Letn^^^g" 
be the projections associated to this .splitting and Ug''^'^{t) = [/e(i)|f|;^^^ • Moreover 
there exist /3i,2,3 > and Ch > independent of 9 satisfying P^ < Pi, P^ < P2 and 
such that the splitting is characterized by the following rate conditions: 

\\U^{t)U^{ry^\\p < Che~^'^'-^\ t>T>0, 
\\U;^me{r)-'\\p < Che^^^'--\ t < r < 0, 
\\U^m^,{T)-'\l<Che^^^'-^\, t,reR 
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and 

(175) \\n*^il'--u''^^^^\\,<cd-^E\i, 

(176) \f3: -(3^\<CS-^\\E\\p, z= 1,2,3, 

(177) C*^^Ch. 

Proof. We make the same identifications as in the proof of CoroUarv l8.21l Consider 
the auxiliary hnear equation 

(178) A{t)^Al{e + ojt)A{t) 
with 

'AlHe) 

22/ 




Alie)=\ Af(0) 

V 

Clearly U;{t) = {tj^{t), U^{t), U^{t)) is a solution of (fTTS]) . By hypothesis (6*) - 
j4a(0)|Ip is small. Then the application of Proposition l8.19l gives the results. □ 

Remark 8.23. We can give an alternative proof to Proposition 18.191 parallel to 
the one for maps. We just sketch it for the stable bundle in the following. Recall 
that we have the invariance condition for all times i > 

The graph condition then writes for all times t > 

(m?0)) eGraph(M(0 + u;t)). 

We now consider the time-one map Ui = Ug{l). The graph condition leads to 
a functional equation which is solved by a fixed point argument. To propagate the 
result to any time and get the estimates, one just has to use the co-cycle property 
as stated in Proposition 18. 71 

The other non-degeneracy conditions can be checked in exactly the same way 
(as in the previous section) and we do not repeat the arguments. 

Lemma 8.24. If \\Em~i\\p^_-i^ is small enough, then 

• If DK^_iDKm-i is invertible with inverse N^-i then DK^DKm is in- 
vertible with inverse Nm and we have 

\\N^\\p„. < \\Nra-l\\p,„_, + Cm-lK^6;J^';+'^\\Em-l\\p,„_,. 

• //avg (^m-i) is non-singular then also avg (Sm) is and we have the estimate 

i(avg(^„))-i| < i(avg(5„_i))-'i +c:„-i«'^;/-r'^iii?™-iiip,„-.- 

The last lemma is devoted to the proof of the cohomology obstruction under the 
iterative step. 

Lemma 8.25. Assume \\Em-i\\p,„_i is small enough. If X\^_-^ spans the coho- 
mology of Km^iiT^) at Am-i, then Xx^^ spans the cohomology of KmiT^) at A„i. 
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Proof. We have by assumption that the map 



is an isomorphism. Thanks to the estimates on Am and |Am — Am-i| and the 
continuity of X\ and DX\ with respect to A, we can write 

The previous estimate comes from the identification of the cohomology with the 

_d_ 



)f 

integration over loops of T' and the fact the quantity ■jr-K*ixx^ matrix no- 



tation 

dX 



z,A-,.,-.(«w)(™»). 

map 



This shows the invertibihty of the map 

dA' 

□ 



8.8. Vanishing lemma. This section is devoted to the proof of Theorem 18.101 
First, recaU that the Lie derivative of the 1-form a with respect to a vector-field L 
is given by (Cartan formula) 

Clo; — diLCt + ihda- 
The following result is of general interest and is a vanishing lemma. 

Lemma 8.26. Assume lo e Dh{n, i^) with k > and ly > I — 1, and X\ is a family 
of real analytic symplectic vector- fields. Let K : Dp D T' — > he a solution of 

(179) d^K^XxoK + E, 

for |A — A*| small enough. Assume furthermore that 

(1) Xy* is exact symplectic and for all A G M' and A ^ X* , the vector-fileds X\ 
are symplectic but not exact symplectic. 

(2) For all A G M', X\ can be extended holomorphically to a complex neighbor- 
hood ofK{Dp). 

(3) The family X\ spans the cohomology o/if (T') at X — X* , i.e. the map 

is an isomorphism. 
Then there exists a constant C such that 

\X-X*\ < C\\E\\p. 

Proof. The proof is very similar to the proof of Lemma 14.91 Indeed, if we consider 
vector-fields as "infinitesimal" diffeomorphisms, the present proof can be considered 
as an infinitesimal version of the proof of Lemma [4.91 We define cr^ ^ as in ([M]) . 

dSSl). 

The proof consists of computing 
(181) [ Cx,a 
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in two different ways. First, notice that by Cartan's fomula, we have 



J Koa. i JK 



ix^n. 



Expanding this last expression in terms of A yields 

d 



^x,n^ 1X1^+ U^^d)\x^^.{\-\*) + 0{\\-X*\^). 

■I Koa. a JK^" - UA \ / 



Furthermore, since the vector-field Xx* is exact symplectic, we have Cx^Oi = dW 
and then the first term in the right-hand side vanishes. We are led to 

(182) / ^x.^^ I ^(ix.d\u=x.{\-y) + o{\x~y\^). 

JKoa.g, J Koa. g. "-^^ ' 

On the other hand, using the linearity of the Lie derivative w.r.t. the vector-field 
and equation (|179p . we have 

Cx^a^ / ix^^= / ia^^ + R, 

where R is such that \\R\\p < C\\E\\. 

Furthermore, by the change of variables formula and the exact symplecticness of 
the manifold we have 



i^dK*a. 



Since to is constant, the exterior differentiation commutes with the contraction 
operator and one gets for all 1 < i < ^ 



io^Q — / di^K*a = 0, 

Koa^ "'i 8; 

yielding 

(183) / Cx^a^R. 

We note that z-component of the map (|180p is the integral of 

over the i — th generator of the torus. Then, from (jl82p - ()183p and the implicit 
function theorem, we get the desired result if |A — A*| is small. □ 

We are now in position to prove Theorem 18. 101 



Proof of Theorem \8.10[ We have to introduce a family of vector-fields X\ satisfying 
the non-degeneracy condition (|180p in Lemma 18.261 Since is non-degenerate, 
given a family of closed 1-forms ax such that ctq — and an exact symplectic 
vector-field X there exists a family of symplectic vector-fields Xx such that 
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1) Xo = X. 

2) ix^^ = o-A- 

Condition 2) implies that Xx is indeed symplectic: 

Cxx^ — dix^^ + ixx^^ — d<y\ = 0. 

If we choose ax such that the cohomology class [cta] ^ for A 7^ then Xx will not 
be exact symplectic for A 7^ 0. Indeed, this follows form the calculation 

Lxx^ — dixx'^ + ixxd-a — dix^o^ + ix^^ = dWx + cr^. 

To choose ax consider the torus K{T^). We take a tubular neighborhood N'^{K{T^)) 
of /\(T'). Since it is contractible to X(T'), then H^{K{T^)) - i/i(iV^(X(T'))). 
Now we consider a basis {5j}i<j<i of H^{K{T^)). We define 

i=l 

Then we have 

J J i=i 

Since K is an embedding, {K*5.j}i<j<.i is a basis of H^(T^) and then the map 
V ^ Dx{K*ixx^)v is invertible. □ 

9. Finite-dimensional Hamiltonian flows 

This section is devoted to the application of our method in the Hamiltonian 
vector-field case. In the same spirit as the previous section, one of the motivations 
is the study of Hamiltonian PDE's. 

The result for Hamiltonian flows is based on the study of the equation 

(184) d^K{e) = J{K{e))VH{K{e)), 

where the function H : M. ^ M. \s the Hamiltonian which is supposed to be real 
analytic. 

Equation (|184p expresses the invariance of the range of K under the Hamiltonian 
vector-field Xh — JVH. We assume that A4 is endowed with fl = dx A dy and 
a = —ydx and hence J is constant. Note that the vector-field JVH is exact 
symplectic. Indeed, by definition, we have 

ijVH^ — -dH. 

Then taking W = —H+ij^nct, we have Cxct — dW. More generally, if we consider 
an exact symplectic vector-field X in the sense of Definition 18. li then there exists 
a function H such that X = JVH. 

Consequently, the Hamiltonian framework fits exactly in the exact case as de- 
scribed in the previous section (due to the lack of cohomology obstruction) . How- 
ever since equations of the type (|184p occur in a lot of physical contexts, our mo- 
tivation to write this section is to provide the formulas showing up for this type of 
system. 

Again, the linearized equation 

(185) ^ = JDVH(K(9 + Lot))A 

dt 
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plays a crucial role. Since A{9) = JD\7H{K{e)) is bounded, equation plS]) admits 
an evolution operator, denoted Ug{t). We have 

j^Ueit)^Ai0 + Lot)Ue{t), 
and C/e(0) = Id. We now have the following definitions. 

Condition 9.1. • Spectral conditions: the evolution operator Ug{t) sat- 

isfies the non-degeneracy Conditions \8.4\ 

• Twist condition: letN{e) = [DK{eyDK{0)]-^ andP{e) = DK{9)N{0). 
The average on T' of the matrix 

S{9) = N{e)DK{9y[A{d)J - JA{9)]DK{9)N{e). 

is non-singular. Here A{d) = JD\7H{K{9)). 

For the sake of completeness, we state a theorem for equation (|184p . It provides 
the existence of invariant tori. 

Remark 9.2. To obtain the expression for S, we used the fact that DK^ JDK = 
and J^^ — — J. 

Theorem 9.3. Let u> satisfy the Diophantine Condition \2.^ Assume the following 
hypotheses 

• The embedding Kq satisfies the non-degeneracy Condition ] 9. 11 

• The map H is real analytic and it can be extended holomorphically to some 
complex neighborhood of the image under Kq of Dp^ : 

Br^{ze C^'^l 39 e {|Im 6*1 < po} s.t. \z - Ko{9)\ < r} , 

for some r > 0. 
Define the error Eq by 

Eo = d^Ko{9) - J\7H{Ko{9)). 

There exists a constant C > depending on I, k, v, \H\(j3[Br)' ll-^-^ollpoi ll-^ollpo^ 
||S'o||po, |(avg (5o))^"'^ I (where Sq and Nq are as in Condition \9. ll revlacina K by Kq) 
and the norms of the projections I|nxo'(e) Upo such that, if Eq satisfies the estimates 

Ck^S-^^'WEqWp, < 1 

and 

CK'S~'''\\Eo\\p„<r, 

where < S < min(l,po/12) is fixed, then there exists an embedding Koo and 
Aoo € such that (Aoo, Koo) G ND{poo ■= Po — 6(5) and 

(186) d^Koo{9) = JVH{Koo{9)). 

Furthermore, we have the estimate 

~ K4p^ <CkH-^^\\E4p,. 

Remark 9.4. One could also formulate a local uniqueness result in the case of 
vector-fields. 
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